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ABSTRACT 

We present a detailed analysis of the X-ray properties of the cooling flows in a sample of 
nearby, X-ray bright clusters of galaxies using high-quality ASCA spectra and ROSAT 
X-ray images. We demonstrate the need for multiphase models to consistently explain 
the spectral and imaging X-ray data for the clusters. The mass deposition rates of the 
cooling flows, independently determined from the ASCA spectra and ROSAT images, 
exhibit reasonable agreement. We conflrm the presence of intrinsic X-ray absorption in 
the clusters using a variety of spectral models. We also report detections of extended 
100/xm infrared emission, spatially coincident with the cooling flows, in several of the 
systems studied. The observed infrared fluxes and flux limits are in good agreement 
with the predicted values due to reprocessed X-ray emission from the cooling flows. 
We present precise measurements of the abundances of iron, magnesium, silicon and 
sulphur in the central regions of the Virgo and Centaurus clusters. Our results firmly 
favour models in which a high mass fraction (70 — 80 per cent) of the iron in the 
X-ray gas in these regions is due to Type la supernovae. Finally, we present a series 
of methods which may be used to measure the ages of cooling flows from the X-ray 
data. The results for the present sample of clusters indicate ages of between 2.5 and 
7 Gyr. If the ages of cooling flows are primarily set by subcluster merger events, then 
our results suggest that in the largest clusters, mergers with subclusters with masses 
of ^ 30 per cent of the flnal cluster mass are likely to disrupt cooling flows. 

Key words: galaxies: clusters: general - cooling flows - intergalactic medium - 
X-rays: galaxies 



1 INTRODUCTION 

In the central regions of most clusters of galaxies, the cooling 
time of the hot intracluster medium is significantly shorter 
than a Hubble time (Edge, Stewart & Fabian 1992; White, 
Jones & Forman 1997; Peres et al. 1998). This is normally 
taken to indicate the presence of a cooling fiow, in which 
the hot gas cools and is deposited throughout the cluster 
core (see Fabian 1994 for a review). Cooling flows are a phe- 
nomenon most easily studied at X-ray wavelengths, where 
the bulk of the cooling occurs. However, large Faraday rota- 
tion measures, consistent with cooling-flow models, are also 
detected in radio observations of cluster cores (Taylor, Bar- 
ton & Ge 1994; Taylor, Allen & Fabian 1999) and strong 
optical emission-line nebulosities and UV/blue emission, as- 
sociated with young, massive stars, are commonly observed 
in and around the dominant galaxies at the centres of cluster 
cooling-fiows {e.g. Johnstone, Fabian & Nulsen 1987; Mc- 



Namara & O'Connell 1989; Allen 1995; Cardiel et al. 1995, 
1998; Crawford et al. 1999). The mass in young stars within 
the central 5 — 20kpc regions of cooling-flow clusters typi- 
cally accounts for < 10 per cent of the total mass inferred to 
be deposited by the cooling fiows within their cooling radii 
(fcooi ~ 50 — 200kpc) although can account for a significant 
fraction of the material deposited within the innermost part. 
At present, however, the fate of the cooling gas at larger radii 
remains uncertain (Daines, Fabian & Thomas 1994; Fabian, 
Johnstone & Daines 1994a; O'Dea et al. 1994; Voit & Don- 
ahue 1995; O'Dea & Baum 1996; Henkel & Wiklind 1998; 
Allen 2000; Ferland, Fabian & Johnstone 2000). 

The X-ray surface brightness profiles of cooling-fiow 
clusters require that the deposition of cooled matter is typ- 
ically distributed throughout the cooling radii with M ^ r. 
This requires that cooling flows are multiphase (Nulsen 
1986; Thomas, Fabian & Nulsen 1987). Spatially- resolved 
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X-ray spectroscopy, carried out with the Position Sensi- 
tive Proportional Counter (PSPC) on ROSAT, has shown 
that the mean emission-weighted temperatures of cooling- 
flow clusters decrease towards their centres (e.g. Allen & 
Fabian 1994; Nulsen & Bohringer 1995). More detailed mul- 
tiphase spectral studies with ASCA have also revealed the 
presence of at least two gas phases in the central regions 
of cooling-flow clusters, with a ratio of emission measures 
in good agreement with the predictions from cooling-flow 
models {e.g. Fukazawa et al. 1994; Fabian et al. 1994b; Allen 
2000). 

Early studies based on ASCA spectra for small samples 
of bright, nearby cooling flows {e.g. Fabian et al. 1994b) 
were complicated by uncertainties associated with the cal- 
culation of the Fe-L emission-line complex used in plasma 
emission codes. These calculations have since been signifi- 
cantly improved in the MEKAL code of Kaastra & Mewe 
(1993; incorporating the Fe L calculations of Liedhal, Os- 
terheld & Goldstein 1995, although residual uncertainties in 
the calculations of other potentially important elements re- 
main) . Further complications in the modelling of the X-ray 
spectra arise from abundance gradients {e.g. Fukazawa et al. 
1994; Matsumoto et al. 1996; Allen & Fabian 1998), the ef- 
fects of intrinsic absorption {e.g. White et al. 1991; Allen & 
Fabian 1997; Wise & Sarazin 2000) and temperature varia- 
tions within the flows, which depend upon the details of the 
local gravitational potentials. 

In this paper we study ASCA spectra and ROSAT X- 
ray imaging data for seven bright, nearby cooling flows; the 
Perseus Cluster (Abell 426), Abell 496, the Virgo Cluster, 
the Centaurus Cluster (Abell 3526), Abell 1795, 2199 and 
2597. We also include one non-cooling flow cluster, the Coma 
Cluster (Abell 1656), for comparison purposes. Our primary 
goal is to determine the properties of the cooling flows, us- 
ing the improved MEKAL code in both our spectral and 
imaging analyses. We examine regions of the clusters large 
enough to encompass the entire cooling flows, thereby min- 
imizing systematic problems associated with modelling the 
complex point spread functions of the ASCA mirrors and 
the underlying spatial distributions of the X-ray emitting 
and absorbing gas. 

Section 2 of this paper describes the observations and 
data reduction. Sections 3 and 4 discuss the method and ba- 
sic results from the spectral analysis. Section 5 examines the 
enrichment history of the ICM in the central regions of the 
Virgo and Centaurus clusters. Section 6 describes the de- 
projection analysis of the imaging data. Section 7 presents 
a comparison of the properties of the cooling flows deter- 
mined from the independent spectral and imaging methods. 
The evidence for intrinsic X-ray absorption and reprocessed 
emission at infrared wavelengths are explored. Section 8 ex- 
amines the constraints that may be placed on the ages of 
cooling flows from the X-ray data. Section 9 discusses the 
possible role of heating process in cooling flows and the likely 
effects of subcluster merger events. The most important 
conclusions from our work are summarized in Section 10. 
Throughout this paper, we assume ^0=50 kms~^ Mpc~^, 
Q. = 1 and A = 0. For the Virgo Cluster, a distance of 18Mpc 
is assumed. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The ASCA observations 

The ASCA (Tanaka, Inoue & Holt 1994) observations were 
made over a three-and-a-half year period between 1993 
May and 1996 December. The ASCA X-ray telescope ar- 
ray (XRT) consists of four nested-foil telescopes, each fo- 
cussed onto one of four detectors; two X-ray CCD cam- 
eras, the Solid-state Imaging Spectrometers (SIS; SO and 
SI) and two Gas scintillation Imaging Spectrometers (CIS; 
G2 and G3). The XRT provides a spatial resolution of ~ 3 
arcmin (Half Power Diameter) in the energy range 0.3 — 12 
keV. The SIS detectors provide excellent spectral resolution 
[AE/E = 0.02(_E/5.9keV)"°-^] over a 22 x 22 arcmin^ fleld 
of view. The CIS detectors provide poorer energy resolution 
[A_B/£ = 0.08(S/5.9keV)"°-^] but cover a larger circular 
field of view of ~ 50 arcmin diameter. 

For our analysis of the ASCA data we have used the 
screened event lists from the revl processing of the data sets 
available on the GSFC ASCA archive (for a detailed descrip- 
tion of the revl processing see the GSFC ASCA Data Re- 
duction Guide, published by GSFC.) The ASCA data were 
reduced using the FTOOLS software (version 3.6) issued 
by GSFC, from within the XSELECT environment (ver- 
sion 1.3). Further data-cleaning procedures as recommended 
in the ASCA Data Reduction Guide, including appropriate 
grade selection, gain corrections and manual screening based 
on the individual instrument light curves, were followed. A 
summary of the ASCA observations, including the individ- 
ual instrument exposure times after all screening procedures 
were carried out, is given in Table 1. 

Spectra were extracted from all four ASCA detectors 
(except for Abell 1795 where the SI data were lost due to sat- 
uration problems caused by fiickering and 'hot' pixels in the 
CCDs). The spectra were extracted from circular regions, 
centred on the peaks of the X-ray emission. For the SIS 
data, the radii of the regions used were selected to minimize 
the number of chip boundaries crossed (thereby minimizing 
the systematic uncertainties introduced by such crossings) 
whilst covering as large a region of the clusters as possible. 
Data from the regions between the chips were masked out 
and excluded. The final extraction radii for the SIS data 
are summarized in Table 2. The Table also notes the chip 
modes used for the observations (1,2 or 4 chip mode) and the 
number of chips from which the extracted data were drawn. 
For the CIS data, a fixed extraction radius of 6 arcmin was 
adopted. 

Background subtraction was carried out using the 
'blank sky' observations of high Galactic latitude fields com- 
piled during the performance verification stage of the ASCA 
mission. The background data were screened and grade se- 
lected in the same manner as the target observations and 
the background spectra were extracted from the same re- 
gions of the detectors as the cluster spectra. (We assume 
that the errors associated with the background subtraction 
are statistical in origin.) For the SIS data, response matri- 
ces were generated using the FTOOLS SISRMG software. 
Where the spectra covered more than one chip, response 
matrices were created for each chip and combined to form a 
counts-weighted mean matrix. For the CIS analysis, the re- 
sponse matrices issued by GSFC on 1995 March 6 were used. 
For both the SIS and CIS data, auxiliary response files were 
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generated with the ASCAARF software, with effective area 
calculations appropriate for extended sources. 



2.2 The ROSAT Observations 

The ROSAT observations were carried out between 1991 
February and 1994 August. For the more distant systems 
at redshifts z > 0.03, the observations were made with the 
High Resolution Imager (HRf), which provides a ~ 5 arcsec 
(FWHM) X-ray imaging facility covering a ~ 40x40 arcmin^ 
field of view (David et al. 1996). For the clusters at lower 
redshifts, the PSPC was used, which has a lower angular 
resolution (~ 25 arcsec FWHM) but provides a larger (2x2 
degree^) field of view and enhanced sensitivity to extended, 
low-surface brightness X-ray emission. Using this combina- 
tion of detectors both good spatial resolution and a field of 
view extending beyond the cooling radii in all clusters was 
obtained. 

The reduction of the data was carried out using the 
Starlink ASTERfX package. For the HRI data. X-ray im- 
ages were extracted on a 2 x 2 arcsec^ pixel scale. For the 
PSPC observations, only data in the 0.4 — 2.0 keV band 
were used, and the images were extracted with a pixel scale 
of 15 X 15 arcsec^. Where more than one observation of a 
source was made, a mosaic was constructed from the indi- 
vidual observations. In each case the data were cleaned and 
corrected for telescope vignetting and accurate centres were 
determined for the X-ray emission from the clusters. A sum- 
mary of the ROSAT observations and the coordinates of the 
X-ray centroids is supplied in Table 3. 



3 SPECTRAL ANALYSIS OF THE ASCA DATA 

3.1 The basic spectral models 

The modelling of the X-ray spectra has been carried out us- 
ing the XSPEC spectral fitting package (version 9.0; Arnaud 
1996). For the SIS data, only counts in pulse height analyser 
(PHA) channels corresponding to energies between 0.6 and 
10.0 keV were included in the analysis (the energy range 
over which the calibration of the SIS instruments is best- 
understood). For the GIS data, only counts in the energy 
range 1.0 — 10.0 keV were used. The spectra were grouped 
before fitting to ensure a minimum of 20 counts per PHA 
channel, allowing x^ statistics to be used. 

The spectra have been modelled using the plasma codes 
of Kaastra &: Mewe (1993; incorporating the Fe L calcula- 
tions of Liedhal, Osterheld & Goldstein 1995) and the photo- 
electric absorption models of Balucinska-Church & McCam- 
mon (1992). The data from all four detectors were included 
in the analysis, with the fit parameters linked to take the 
same values across the data sets. The exceptions to this were 
the emission measures of the ambient cluster gas in the four 
detectors which, due to the different extraction radii used 
and residual uncertainties in the flux calibration of the in- 
struments, were maintained as independent fit parameters. 

The spectra were first examined with a series of four, 
basic spectral models: model A consisted of an isothermal 
plasma in collisional equilibrium at the optically-determined 
redshift for the cluster, and absorbed by the nominal Galac- 
tic column density (Dickey & Lockman 1990). The free 



parameters in this model were the temperature (kT) and 
metallicity [Z) of the plasma (measured relative to the so- 
lar photospheric values of Anders & Grevesse 1989, with the 
various elements assumed to be present in their solar ratios) 
and the emission measures in each of the four detectors. The 
second model, model B, was the same as model A but with 
the absorbing column density (A'h) also included as a free 
parameter in the fits. The third model, model C, included 
an additional component explicitly accounting for the emis- 
sion from the cooling flows in the clusters. The material in 
the cooling flows was assumed to cool at constant pressure 
from the ambient cluster temperature, following the pre- 
scription of Johnstone et al. (1992). The normalization of 
the cooling-flow component was parameterized in terms of 
a mass deposition rate (Ms), which was a free parameter in 
the flts (linked to take the same value in all four detectors) . 
The metallicity of the cooling gas was assumed to be equal 
to that of the ambient cluster gas. The emission from the 
cooling flow was also assumed to be absorbed by an intrin- 
sic column density, AA'^e, of cold gas, which was a further 
free parameter in the fits. The abundances of metals in the 
absorbing material were fixed to their solar values (Anders 
& Grevesse 1989). 

The fourth spectral model, model D, was similar to 
model C except that the constant-pressure cooling fiow com- 
ponent was replaced with a second isothermal emission com- 
ponent. Model D constitutes a more general model which 
should normally provide an adequate description of the more 
specific cooling-flow models at the spectral resolution and 
signal-to-noise ratios typical of ASCA observations (see dis- 
cussion in Section 3.4). The temperature and normalization 
of the second, cooler emission component were included as 
free parameters in the fits (although the normalization of 
the second emission component was linked to take the same 
value in all four detectors) . The second emission component 
was again assumed to be intrinsically absorbed by a col- 
umn density, AA*'h, of cold gas, which was a further free 
parameter in the fits. The metallicities of the two emission 
components were assumed to be equal. 

Fig. 1 shows the ASCA data and best-fitting spectral 
models for four of the clusters in the sample; Abell 426, 
the Virgo cluster, the Centaurus cluster and Abell 2199. 
For clarity, only the results for the SI and G2 detectors are 
shown. Table 4 summarizes the fit results for the clusters, 
using the four basic spectral models. We note that the SO and 
G3 data for the Centaurus Cluster and the SO and SI data 
for Abell 2597 were found to exhibit small (~ 1 per cent) 
gain mis-calibrations, for which corrections were applied. 

3.2 Isothermal cooling flow models 

The constant pressure cooling flow model (model C) pro- 
vides a crude approximation to the X-ray emission spectrum 
from a multiphase cooling flow. Recently, Nulsen (1998) has 
presented a more sophisticated, self-similar model for the 
emission from a cooling flow in which the mean gas tem- 
perature remains constant with radius i.e. an 'isothermal' 
cooling flow. Although the isothermal cooling-flow model is 
rather speciflc, the assumption of an approximately constant 
mass-weighted temperature over the bulk of the radius in a 
cluster cooling flow is consistent with current constraints on 
the distribution of gravitating matter in the core regions of 
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massive clusters from X-ray and gravitational lensing studies 
(Allen 1998) and with previous, spatially-resolved spectro- 
scopic X-ray studies of nearby clusters cooling flows with 
ASCA {e.g. Fukazawa et al. 1994; Ohashi et al. 1997; Ikebe 
et al. 1999). The assumption of an approximately isother- 
mal mean gas temperature in cooling flows also leads to a 
more plausible implied distribution of initial density inho- 
mogeneities in the cluster gas than would be the case if the 
temperature profiles dropped significantly within the core 
regions (Thomas, Fabian & Nulsen 1987). 

We have incorporated the Nulsen (1998) isothermal 
cooling flow model into the XSPEC code and applied it 
to the ASCA data. As with the constant pressure cooling 
flow model (model C), the isothermal cooling-flow model 
introduces only two extra free parameters into the fits over 
and above those present in the single-phase spectral model 
(model B); the mass deposition rate (Ms) and the column 
density of intrinsic X-ray absorbing material (AA^h). Note, 
however, that whereas in the case of the constant pressure 
cooling flow model the tabulated temperatures are the up- 
per temperatures from which the gas cools, for the isother- 
mal cooling flow models the quoted temperatures are the 
mean temperatures (T from the work of Nulsen 1998) in the 
cooling flows, which are directly related to the gravitational 
potentials. (For a comparison of the emission spectra from 
isothermal and constant pressure cooling flow models see 
Nulsen 1998). 

For our initial analysis, we assumed that the slope of the 
mass deposition profile, rj = 1 (where M oc r**). The mean 
temperature of the gas in the isothermal cooling flows was 
linked to be equal to that of the ambient gas at larger radii 
in the clusters (which, as with spectral models A-D, was a 
free fit parameter) . The results determined with the isother- 
mal cooling fiow model with rj = 1 (hereafter spectral model 
E) are listed in Table 4. In most cases, the isothermal cool- 
ing flow model provides at least as good a fit to the ASCA 
spectra as the constant pressure model (model C) , although 
for the Virgo cluster spectral model C provides a better fit. 
We note that for the Virgo cluster, the central ~ 5 arcmin 
radius (SO) aperture analysed corresponds to a spatial scale 
of only ~ 25 kpc. For this system, it seems plausible that 
the gravitational potential may flatten significantly in the 
region studied, and that the approximation of isothermal- 
ity in the mass-weighted X-ray temperature profile may no 
longer hold. 

We have also investigated the effects on the spectral fits 
of using other values of ?;, in the range 0.75 — 2.5. These re- 
sults are shown in Table 5. As discussed in Section 8.3, a 
value of 77 ^ 1.5 provides a good match to the image depro- 
jection analysis results for a number of the clusters studied, 
with statistically acceptable values spanning the range 0.75 
to 2.5. The effect of increasing the value of 77 in the isother- 
mal cooling flow models is always to increase the value of 
the mass deposition rate and to increase the intrinsic absorp- 
tion required on that component. The Virgo cluster shows 
the largest increase in mass deposition rate; a factor of two 
when increasing r) from 0.75 to 2.5. Other clusters show more 
moderate increases of 20 — 70 per cent. In general, the value 
of x^ also decreases as r] is increased. (The decrease in some 
cases is statistically significant when increasing rj from 0.75 
to 2.5.) Abell 426 is the only cluster for which a minimum 
in x^ is found. In most cases the reduction in y^ continues 



to implausibly high values of rj. We do not find a signifi- 
cant improvement in x^ when increasing rj from the value 
rj — 1.0 used in model E to the value of rj = 1.5, the most 
likely value for several of the clusters determined from the 
imaging data. 

In what follows, we will adopt the results determined 
with spectral model E (77 = 1) as our 'best' spectral results 
on the cooling fiows, except for the Virgo Cluster for which 
the results determined with the constant pressure model 
(model C) are preferred. 

3.3 The requirement for multiphase models and 
excess absorption 

The results determined from the basic spectral analysis, 
summarized in Table 4, show that (with the exception of 
Abell 1795) the fits with the single-temperature models (A 
and B) are significantly improved (i.e. a significant reduc- 
tion in x^ is obtained) when the absorbing column density 
acting along the line of sight is included as a free param- 
eter in the fits. (For guidance, a reduction of Ax^ ~ 7, in 
a fit with u ~ 1000 degrees of freedom and a reduced chi- 
squared value xi "^ 1.00, is significant at approximately the 
99 per cent confidence level.) A further highly significant 
reduction in x^ is obtained for all clusters when including 
either an absorbed cooling flow spectrum (models C, E) or 
an absorbed, cooler isothermal emission component (model 
D) in the modelling. (For guidance, improvements in x^ of 
~ 10 and 15, with the introduction of 2 (models C, E) and 
3 (model D) extra fit parameters, are significant at approxi- 
mately the 99 per cent confidence level in a fit with i' ~ 1000 
and xl ~ 1-00.) 

The requirement for the introduction of cooler omis- 
sion and absorption components is also illustrated in Fig. 
2, where we show the residuals (the data divided by the 
model predictions, in units of x) from fits to the individ- 
ual SIS spectra at energies above 3 keV, using a simple, 
single temperature model with Galactic absorption (model 
A). The best-fitting models have then been extrapolated to 
cover the full 0.6 — 10.0 keV band of the SIS detectors. In 
all cases, except the non cooling-flow Coma cluster, a clear 
excess in the residuals is observed at energies between ~ 0.8 
and 3.0 keV. This region of the spectrum is dominated by 
the Fe-L, Mg, Si and S line complexes, and the detection of 
positive residuals due to these features provides a clear in- 
dication of the presence of gas significantly cooler than the 
ambient cluster temperature. In addition, the residuals for 
a number of the clusters exhibit a clear deficit at low ener- 
gies {E < 0.8 keV) providing evidence for excess absorption, 
over and above that due to our own Galaxy, associated with 
the central regions of the clusters. 

The temperatures and metallicities determined from the 
fits to the SIS data in the 3—10 keV energy range are sum- 
marized in Table 6. The results are in good agreement with 
those determined from the more-sophisticated, multiphase 
analyses of the full ASCA data sets, using spectral models 
C, D and E. 

The residual diagram for Abell 426, determined from 
the 3 — 10 keV SIS data using the nominal (Dickey & Lock- 
man 1990) value for the Galactic column density (A''h = 
1.49 X lO'^^ atom cm^^), shows a strong excess of emission 
below an energy, E ~ 0.8 kcV (see the dotted points in Fig. 
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2). We suspect this excess to be due to an over-estimation 
of the Galactic column density in the region covered by the 
ASCA spectra, which is significantly smaller than the resolu- 
tion of the HI studies (Dickey & Lockman 1990). In order for 
the residual diagram for Abell 426 to appear like those for 
the other clusters included in this study, a Galactic column 
density of ~ 1.0 x 10^^ atom cm^^ is required (solid points 
in Fig. 2). The ROSAT PSPC X-ray colour ratio study of 
the Perseus cluster by Allen & Fabian (1997) also suggests 
a Galactic column density in the region of the cluster core 
of ~ 1.1 X 10^^ cm"^ . The lOO/im IRAS map of the Perseus 
cluster is complex (the system lies at a Galactic latitude of 
6 = — 13 degrees) and allows for a decrease in the Galactic 
column density towards the cluster centre. The result on the 
reduced Galactic column density for Abell 426 is also sup- 
ported by the detailed modelling of the full ASCA data sets 
with spectral models C, D and E, which determine values 
in the range 0.8 — 1.2 x 10^^ atom cm~^. The results on the 
Galactic column densities for the other clusters are in rea- 
sonable agreement with the values determined from the HI 
studies (Dickey & Lockman 1990). 

3.4 A note on two-temperature models 

The results listed in Table 4 show that the two-temperature 
spectral model (model D) generally provides a slightly better 
description of the ASCA spectra than the cooling-flow mod- 
els (C, E). This has sometimes been taken to indicate that 
the cluster gas is distinctly two-phase, with the cooler phase 
being due to the central galaxy {e.g. Makishima 1997; Ikebe 
et al. 1999). However, the longevity of such a gas phase im- 
plied by the common occurrence of cooling flows would then 
require the presence of a heat source to carefully balance the 
radiative cooling. The very high cooling rates found in some 
distant, luminous clusters {e.g. Allen 2000) are also difficult 
to explain with such models. 

The fact that the two-temperature models often pro- 
vide a better description of the ASCA data is not, a pri- 
ori, surprising. Model D provides a more general fitting pa- 
rameterization for multi-temperature plasmas observed with 
ASCA, with an extra degree of freedom over the cooling flow 
models. Simulated cooling-flow spectra, constructed at the 
spectral resolution and with count statistics comparable to 
those of typical ASCA observations, are well-described by 
two-temperature models {e.g. Allen 2000; see also Buote, 
Canizares & Fabian 1999). Where spectral model D pro- 
vides a significantly better fit to the observations than mod- 
els C and E, this need only indicate that the simple cooling 
flow models over-simplify the true emission spectra from the 
central regions of the clusters. In particular, given that the 
systems studied here are the nearest, X-ray brightest cluster 
cooling flows, in which strong radial variations in metallic- 
ity and intrinsic absorption are known exist (Fukazawa et al. 
1994, Matsumoto et al. 1996, Allen & Fabian 1997) and in 
which clear deviations from solar abundance ratios are ob- 
served (Section 5), the fact that the basic cooling-flow mod- 
els over-simplify the observed situation is to be expected. 

We note that for the very best ASCA data sets, such 
as the observation of the Centaurus Cluster presented here, 
the data are of sufficient quality to demonstrate the presence 
of more than two thermal components. Section 4 discusses 
the evidence for additional power-law components in the fits 



with the two-temperature plasma models. For the Centau- 
rus Cluster, the fit with the two-temperature plus power-law 
model (also incorporating variable element abundance ra- 
tios; see Section 5) is further improved by the introduction 
of a third plasma component. Following Allen, Di Matteo 
& Fabian (2000) and using only the data from the SI and 
G2 detectors for the fits incorporating the power-law com- 
ponents (which are the best calibrated detectors for that 
particular observation, the SO and G3 data sets exhibiting 
noticeable gain offsets; Section 3.1), we find that the intro- 
duction of the third plasma component (which is assumed 
to be absorbed by the same intrinsic column density as the 
power-law and the cooler of the two initial plasma compo- 
nents) leads to a drop in x^ of Ax^ = 37 for the introduction 
of 2 additional degrees of freedom. The best fit temperatures 
in this case are 3.25, 1.51 and 0.64 keV, respectively. 

3.5 Comparison with previous results 

A number of previous works have also presented results 
based on ASCA observations of the central regions of the 
clusters studied here. Fabian et al. (1994b) report results 
from an early analysis of ASCA data for the innermost re- 
gions of the Perseus (Abell 426) , Centaurus and Coma clus- 
ters and Abell 1795. The regions studied by Fabian et al. 
(1994b) were smaller than those studied here and did not 
include the the entire cooling flows. For this reason the mass 
deposition rates measured by Fabian et al. (1994b) are ap- 
proximately a factor two smaller than those reported here. 

Fukazawa et al. (1994) present an analysis of a previous, 
shorter ASCA observation of central regions of the Centau- 
rus Cluster (the same data studied by Fabian et al. 1994b). 
These authors reported the discovery of a strong metallic- 
ity gradient in the cluster, with a central metallicity and 
temperature in good agreement with the results presented 
here. Ikebe et al. (1999) present results from a detailed anal- 
ysis of the longer, more recent observation of the Centau- 
rus Cluster also studied in this paper. Their results from 
fits with two-temperature models incorporating variable el- 
ement abundance ratios are in good general agreement with 
those reported here (Section 5). 

Matsumoto et al. (1996) present an analysis of ASCA 
observations of the central regions of the Virgo Cluster. 
Their results for the innermost regions of the cluster, de- 
termined using a two temperature model with variable el- 
ement abundances, are in reasonable agreement with those 
reported here, although small systematic differences exist. 
These differences are presumably due to the different plasma 
codes used in the two studies. The MEKAL code (Kaastra 
& Mewe 1993; Liedhal, Osterheld & Goldstein 1995) em- 
ployed in this paper contains significant improvements over 
the Raymond & Smith (1977) code used in the earlier Mat- 
sumoto et al. (1996) work. Our results on the temperature 
and abundance structure in the central regions of the Virgo 
and Centaurus clusters are in good agreement with the re- 
cent results of Buote et al. (1999). 

Xu et al. (1998) report results from a study of the cen- 
tral regions of AbeU 1795 with ASCA. Their results, ob- 
tained using a simple, two-temperature model are in good 
agreement with those reported here from our analysis with 
spectral model D. The temperature in the central regions 
of the Coma cluster determined with spectral model B is 
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in good agreement with the previous result of Honda et al. 
(1996). 

Markevitch et al. (1998, 1999) present detailed results 
on the temperature structure in Abell 496, 1795, 2199 and 
2597. These authors determine mean temperatures, cor- 
rected for the effects of cooling flows, for these clusters of 
4.7±0.2, 7.8±1.0, 4.8±0.2 and 4.4l°-j keV, respectively. In 
general, the mean cluster temperatures measured by Marke- 
vitch et al. (1998, 1999) are slightly hotter than the results 
for the central regions of the clusters reported here, using 
spectral models C and E (the cooling- flow models). 



4 THE EVIDENCE FOR POWER-LAW 
COMPONENTS 

Allen et al. (2000) and Di Matteo et al. (2000) have pre- 
viously reported the detection of hard, power-law emission 
components in the ASCA data for the Virgo and Centaurus 
clusters. The presence of these components is clearly evident 
in the residual diagrams shown in Fig. 2. We have extended 
this analysis to the other higher luminosity systems reported 
here. Starting in each case with the best-fltting two temper- 
ature model, we have examined the improvements to the 
fits obtained by introducing an additional power-law emis- 
sion component.n The photon index and normalization of 
the power-law components were included as free fit param- 
eters. (Note that for the Virgo and Centaurus clusters, the 
best fitting two-temperature models incorporating variable 
element abundances were used as the starting points for the 
analysis; see Section 5). The power-law components were 
assumed to be absorbed by the same intrinsic column den- 
sity as the cooler thermal emission components.r|The results 
from the fits are summarized in Table 7. 

The results for Abell 426 and 2597 suggest the presence 
of relatively luminous power-law components (Lx.i-io = 



A.Otli X 10" 



and 2.0^og x 10''"' ergs ^, respectively) 



* We use the two-temperature model (model D) as the start- 
ing point for our analysis of the power-law components since it 
provides the most flexible fitting parameterization for the cluster 
emission, and therefore the most rigorous test for the presence of 
additional power-law components. The two-temperature model 
also allows for the inclusion of variable element abundance ratios 
(using the vmekal model in XSPEC), which is firmly required 
by the data for the Virgo and Centaurus clusters (Section 5). We 
note that the inclusion of a power-law component into the fits with 
the cooling fiow emission models (models C and E) also leads to 
significant improvements in x^ and constraints on the power-law 
components consistent with those listed in Table 7 (although for 
Abell 496 the inferred 2—10 keV flux is approximately twice as 
high and the x^ value +46 worse). For the Virgo and Centaurus 
clusters, Abell 1795 and 2199, the mass deposition rates deter- 
mined with the cooling fiow-plus-power-law models are in good 
agreement with the values listed in Table 4, although for Abell 
426, 496 and 2597, the rates are reduced by a factor > 1.6. 
1 The ASCA spectra cannot easily constrain more complex spec- 
tral models in which the intrinsic column density acting on the 
power-law component is included as a separate free parameter, 
although such models are probably more realistic. The intrinsic 
column densities associated with the power-law components mea- 
sured with such models are, however, likely to be larger than the 
values listed in Table 7). 



with photon indices (F = 2.05^0.05 and l.SOto.tl) con- 
sistent with the typical values determined for Seyfert nu- 
clei {e.g. Nandra et al. 1997). However, the flux associated 
with the power-law component in Abell 426 (i^x,2-io = 
1.80^0 14 X 10"^" ergcm"^ s~^ ) is approximately an order 
of magnitude larger than the nuclear point-source flux in- 
ferred from ROSAT HRI observations of the dominant clus- 
ter galaxy, NGC 1275, made in 1991 February and 1994 
August (using the same spectral model). This implies ei- 
ther the presence of signiflcant intrinsic absorption associ- 
ated with the nucleus of NGC 1275 (see also Section 7.3), 
that the power-law emission originates (at least in part) from 
some other source, or some other deflciency with the spectral 
model. 

The results for Abell 496, 1795 and 2199 are more anaF 
ogous to those those previously reported for the Virgo and 
Centaurus clusters. For Abell 496, we measure a photon in- 
dex, F = 1.44^033 and an intrinsic 1 — 10 keV luminosity 
associated with the power-law component of 3.29^084 x 10*^ 
ergs~^ , which is ~ 80 times larger than the value for M87. 
The results for Abell 1795 and 2199 indicate marginal im- 
provements to the flts with the introduction of hard, power- 
law components, with intrinsic 1 — 10 keV luminosities of ~ 6 
and 3.5 x 10*^ ergs~^ , respectively. In all three cases the 
ASCA data allow for signiflcant intrinsic absorption acting 
on power-law emission, which is consistent with the fact that 
a central point sources appear to contribute only weakly in 
the ROSAT band. 

Comparing the results in Table 7 with those in Table 4, 
we see that the most noticeable effect on the best-fit parame- 
ter values measured with spectral model D, on the introduc- 
tion of the power-law component, is to decrease the temper- 
ature of the hotter thermal component to a value in better 
agreement with the results obtained using the cooling fiow 
models, C and E. For Abell 426, the mean metallicity is also 
increased from ~ 0.5.^0 to 0.8^©, which is consistent with 
the presence of a metallicity gradient in the cluster, albeit 
with the caveats given above. (We note that the relatively 
poor x^ value for Abell 426 is primarily caused by residuals 
in the SO data. If we repeat our analysis of Abell 426 in- 
corporating the power-law component but excluding the SO 
data, we obtain similar best-fit parameters and x^ = 1570 
for 1398 degrees of freedom.) 



4.1 Povirer-law emission from the Coma Cluster 

The spectral fit to the Coma Cluster with spectral model 
D was not improved by the introduction of an additional 
power- law component. We are able to place a 90 per cent 
confidence upper limit on the fiux of any power-law compo- 
nent (with a canonical AGN photon index F = 2.0) incorpo- 
rated into model D of Fx,2-io < 2.2 x 10"^^ ergcm^^ s~^ . 
However, the introduction of a power-law component was 
found to significantly improve the fit obtained with the more 
simple, single-temperature spectral model B {Ax^ = 9.0, 
giving x^ = 1265 for 1225 degrees of freedom). The best-fit 
photon index and normalization for the power-law compo- 
nent, determined using the single-temperature plus power- 
law model, are F = 3.3+?g and Ai = 3.1+i? x 10~^ photon 
keV~^cm~^s~'^. The mean cluster temperature measured 
with this model is 8.5lo'g keV. The 90 per cent confidence 
upper limits on the flux associated with the power-law com- 
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ponent in the 2-10 and 20 - 80 keV bands are 1.9 x 10"" 
and 3.1 x 10"^^ ergcm^^s"^ , respectively. 

Fusco-Femiano et al. (1999) report the detection of 
hard, non-thermal X-ray emission from the the Coma Clus- 
ter using observations made with the Phoswhich Detection 
System (PDS) and High-Pressure Gas Scintillation Propor- 
tional Counter (HPGSPC) on the BeppoSAX satellite. (The 
PDS and HPGSPC detectors have fields of view ~ 1.3 
(hexagonal) and 1.0 degree across, respectively.) From a 
joint fit to the PDS and HPGSPC data, which together 
cover the 4.5 — 80 keV range, using a single-temperature plus 
power-law model, these authors determine a mean temper- 
ature for the ICM of 8.5lg'5 keV, and a best-fitting power- 
law photon index of F = 1.7-3. The 20 - 80 keV flux 
associated with the power-law component (measured from 
the PDS data alone) is ~ 2.2 x 10"" ergcm"^s~^ . These 
parameters are consistent with those determined from the 
ASCA data, using the same simple spectral model. 

Rephaeli, Gruber & Blanco (1999) also report observa- 
tions of the central ~ 1 degree region of the Coma Clus- 
ter using the Proportional Counter Array (PGA) and High- 
Energy X-ray Timing Experiment (HEXTE) on the Rossi 
X-ray Timing Explorer (RXTE). Rephaeh et al. (1999) show 
that the RXTE data in the 3 — 100 keV range can be equally 
well-modelled using either a two-temperature plasma model 
(with temperatures of 8.2^^2 and 4.7^30 keV, respectively; 
with the cooler component accounting for ~ 20 per cent of 
the total flux) or a single-temperature plus power-law model. 
In this second case (favoured by the authors), the RXTE 
data measure a mean cluster temperature of 7.51 ±0.18 keV 
and a photon index and normalization (at 1 keV) associ- 
ated with the power-law component of F = 2.35 ± 0.45 and 
Al = 1.9 ± 0.4 X 10"^ photon keV^^cm^^s"^, respectively. 
This normalization is somewhat higher than that measured 
from the ASCA data, which may indicate that the source 
of the power-law emission in this cluster is more extended 
than the regions included in the ASCA analysis. 



5 THE ABUNDANCES OF ELEMENTS IN 
THE VIRGO AND CENTAURUS CLUSTERS 

5.1 The observed abundances ratios 

The metallicities and abundance ratios of the various ele- 
ments in the ICM place firm constraints on the enrichment 
history of the cluster gas and the relative contributions from 
different types of supernovae. Mushotzky et al. (1996) re- 
ported results from an ASCA study of the outer regions 
of four nearby clusters (Abell 1060, AWM7 and the same 
observations of Abell 496 and 2199 presented here) and con- 
cluded, based on the mean metallicities of O, Ne, Si, S and Fe 
determined from their work, that most of the enrichment in 
the outer regions of clusters is likely to be due to Type H su- 
pernovae. However, these authors warned that a significant 
contribution from Type la supernovae could not be ruled 
out. Subsequent analyses of clusters and groups (Finoguenov 

6 Ponman 1999; Dupke & White 2000a,b) have shown that 
their central regions often exhibit enhanced contributions 
from Type la supernovae. 

The high-quality ASCA spectra for the Gentaurus and 
Virgo clusters analysed here provide an excellent opportu- 



nity to re-examine the issue of individual element abun- 
dances in the ICM in nearby clusters. (The SIS data for 
the Gentaurus and Virgo clusters provide the most signifi- 
cant and reliable constraints of the eight systems examined 
here and we therefore concentrate our discussion on these 
two clusters). Using the SO and SI data for the Virgo clus- 
ter, and the SI data for the Gentaurus Cluster (the data sets 
which are best-calibrated and most sensitive for the present 
task), and starting in each case with the best fitting two- 
temperature model (Model D; which has the abundances 
of all other elements linked to be equal to that of iron), 
we have systematically determined the improvements to the 
fits obtained by allowing the metallicity of each element, in 
turn, to be included as a free parameter in the fit. Having 
identified the element giving the most significant improve- 
ment, the metallicity of that element was maintained as a 
free parameter, and the process repeated to determine the 
element providing the next most significant improvement. 
This was continued until no further significant improvement 
(at > 95 per cent confidence, as determined using an F-test) 
was obtained by allowing other elements to fit freely. 

For the Virgo cluster, the most significant improvement 
to the fit to the joint SO and SI data with spectral model D 
was obtained by including the abundance of Mg as a free fit 
parameter (x^ drops from 757.7 to 570.4). This was followed 
by Si (giving a further drop in x^ to 511.5), Na (x^ ~ 488.8) 
and S (x^ ~ 468.9). At this point including the abundances 
of further elements as free parameters did not significantly 
improve the fit. 

For the Gentaurus cluster, the most significant improve- 
ment in the fit to the SI data was also initially obtained 
by allowing the abundance of Mg to fit freely (x^ drops 
from 478.4 to 409.5). This was followed, in turn, by O 
(X^ = 369.3), Si (x^ = 339.0) and S (x^ = 322.3). 

The results on the abundances of the individual ele- 
ments in the central regions of the Gentaurus and Virgo 
clusters are summarized in Table 8. The two-dimensional 
confidence contours on the abundances of Mg, Si, S and 
Fe in the Virgo Cluster are shown in Fig. 3. The fits with 
the free-fitting element abundances provide highly signif- 
icant improvements in the goodness-of-fit with respect to 
the values determined with the basic two-temperature model 
(model D). However, the best-fitting temperatures {kT and 
kT2), Galactic column densities (A^h) and intrinsic column 
densities (AA'^e) are in good, general agreement with those 
listed in Table 4. The measured Galactic column densities 
are also in good agreement with the values determined from 
HI studies (Dickey & Lockman 1990). We note that the re- 
sults on the element abundances are not significantly modi- 
fied by the inclusion of the power-law components discussed 
in Section 4. 

Our results on the [Si/Fe] and [S/Fe] ratios in the cen- 
tral regions of the Virgo cluster are in good agreement with 
the previous measurements of Matsumoto et al. (1996; al- 
though the absolute element abundances differ, presumably 
due to the different plasma models used in the two studies; 
Section 3.5). Our results on the Si, S and Fe abundances in 
the central regions of the Gentaurus cluster are similar to 
those reported by Ikebe et al. (1999). The measured Si, S 
and Fe abundances in both the Virgo and Gentaurus clus- 
ters are in good agreement with the results of Buote et al. 
(1999). 
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We caution against over-interpretation of the result on 
the low oxygen abundance for the Centaurus Cluster. This 
result is primarily due to features at the lowest energies in 
the SI spectrum, where systematic errors and complexities 
due to absorption in the cluster are at a maximum. Similarly, 
the result on the high Na abundance in the Virgo cluster 
(11 ± 3 times the solar value) is due to a line-like emission 
feature at £ ~ 1.25keV in the spectrum, where residual 
uncertainties in the modelling of the Fe-L emission complex 
remain. 

The element abundances most reliably determined from 
our data are those of Fe, Mg, Si and S. For the Virgo cluster 
we determine abundance ratios (using the solar photospheric 
abundance scale of Anders & Grevesse 1989) of [Mg/Fe]=- 
0.45 (with values in the range -0.67 to -0.28 being consistent 
with the joint 90 per cent confidence limits on the individ- 
ual abundance measurements) , [Si/Fe]=0.18 (0.12-0.24) and 
[S/Fe]=0.10 (0.03-0.17). For the Centaurus cluster we de- 
termine [Mg/Fe]=-0.49 (with a 90 per cent confidence range 
of -0.82 to -0.25), [Si/Fe]=0.14 (0.07-0.22) and [S/Fe]=0.13 
(0.05-0.22). Thus, we find the abundance ratios for the Virgo 
and Centaurus clusters to be both well constrained and in 
excellent agreement with each other. 

Our results on the abundance ratios for [Mg/Fe], [Si/Fe] 
and [S/Fe] in the inner regions of the Centaurus and Virgo 
clusters differ from those of Mushotzky et al. (1996) from 
their study of the outer regions of other nearby clusters 
observed with ASCA. Mushotzky et al. (1996) determined 
mean abundance ratios for [Mg/Fe]=0.07, [Si/Fe]=0.31 and 
[S/Fe]=-0.11 from their sample, using the same solar pho- 
tospheric abundance scale. The differences between these 
results are consistent with the changes in the relative abun- 
dances with radius inferred to be present in other clusters 
{e.g. Dupke & White 2000b). 

5.2 Comparison with supernovae enrichment 
models 

Recently, a number of studies (Ishimaru & Arimoto 1997; 
Gibson, Lowenstein & Mushotzky 1997; Nagataki & Sato 
1998; Dupke & White 2000a) have presented detailed pre- 
dictions for the a— element abundances in the ICM, as a 
function of the fraction of the iron enrichment due to type la 
supernovae, using a variety of theoretical supernovae mod- 
els. Scaling the observed abundance ratios for the central 
regions of the Centaurus and Virgo clusters to the meteoric 
abundance scale of Anders & Grevesse (1989), as used in the 
theoretical calculations, we determine mean, un-weighted 
abundance ratios of [Mg/Fe] ~ -0.63, [Si/Fe] ~ 0.00 and 
[S/Fe]~ 0.02. Thus, the [Si/Fe] and [S/Fe] ratios in the cen- 
tral regions of the Centaurus and Virgo clusters are in good 
agreement with the solar meteoric values. In contrast, the 
observed [Mg/Fe] ratio indicates that magnesium is ~ 4 
times under-abundant with respect to the meteoric value. 

Comparing these results with the predictions of Na- 
gataki & Sato (1998), our observed [Si/Fe] ratio implies a 
mass fraction of the iron enrichment due to type la super- 
novae, MpcSNia/Mpctotai, in the range Mpe.sNia/Mpctotai ~ 
0.6 — 0.9 (where the limits cover the full range of super- 
novae models studied by those authors). A similar com- 
parison of the observed [Si/Fe] result with the range of 
models discussed by Gibson et al. (1997) also requires 



Mpe.sNia/Mpe.totai ~ 0.6 - 0.8. The observed, mean [Mg/Fe] 
ratio favours a mass fraction due to type la supernovae of 
Mpe.sNia/Mpe,totai > 0.7 (Glbsou et al. 1997). 

A comparison of the observed [S/Fe] ratio ([S/Fe]~ 
0.02) with the theoretical models of Nagataki & Sato (1998) 
requires Mpe,sNia/Mpe, total = 0.4 — 0.8, for spherical type 
II supernovae explosions. However, the lower sulphur yields 
associated with aspherical type II explosion models imply 
somewhat lower Mpc,sNia/Mpe,totai values, depending on the 
degree of asphericity and the precise SNIa enrichment model 
used. In particular, the observed [S/Fe] ratio appears incon- 
sistent with enrichment models incorporating both a high 
degree of asphericity in the SNell explosions and simple 
deflagration of the type la supernovae (Nagataki & Sato 
1998). Comparing the observed [S/Fe] value with Gibson et 
al. (1997) indicates Mpc,sNia/Mpe. total ^ 0.7. In contrast to 
previous results (Mushotzky et al. 1996; Gibson et al. 1997), 
our observations do not indicate a severe deficiency of sul- 
phur with respect to the other a-elements, when compared 
with the theoretical predictions (at least for models incor- 
porating relatively low degree of asphericity in the type-II 
supernovae) . 

In summary, the observed a— element abundance ra- 
tios in the central regions of the Centaurus and Virgo 
clusters favour a high mass fraction of the observed iron 
abundance being due to type la supernovae, with a value 
Mpe,sNia/Mpe, total ~ 0.7 — 0.8 providing consistent results. 



5.3 Type la supernovae and the formation of 
metallicity gradients 

Spatially resolved spectroscopy of the Centaurus and Virgo 
clusters has revealed the presence of strong metallicity gra- 
dients in the central regions of these systems (Fukazawa et 
al. 1994, Matsumoto et al. 1996). Similar gradients have also 
been reported in a number of other nearby cooling flows e.g. 
AbeU 496 (Hatsukade et al. 1997) and AWM7 (Ezawa et al. 
1997). Allen & Fabian (1998) showed that the mean emis- 
sion weighted metallicity of luminous cooling-flow clusters is 
approximately twice that of non cooling-flow systems, which 
is most plausibly understood as being due to the presence of 
metallicity gradients in the cooling flow systems. In contrast, 
ASCA observations of nearby, non-cooling flow clusters and 
clusters with relatively small cooling flows e.g. the Coma 
Cluster (Watanabe et al. 1997), AbeU 1060 (Tamura et al. 
1996) and the Ophiuchus Cluster (Matsuzawa et al. 1996), 
show little or no metallicity gradients. 

The results on the individual element abundances in 
the central regions the Centaurus and Virgo Clusters shows 
that, where the metallicity gradients exist, the enrichment of 
the ICM is primarily due to type-la supernovae. The mech- 
anism by which metallicity gradients form in cooling-flow 
clusters remains unclear, but may be due to one or more 
of galaxy stripping, continuous enrichment by supernovae 
within the central galaxy, or a Type la supernova wind from 
the central galaxy (Dupke & White 2000b). The absence of 
metallicity gradients in non-cooling flow systems appears re- 
lated to recent or ongoing subcluster merger events in these 
systems, which will disrupt and mix the central high metal- 
licity gas with the surrounding lower-metallicity material 
(Allen & Fabian 1998). The presence of a metallicity gra- 
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dient then appears a standard property of regular, relaxed 
clusters which also, invariably, contain cooling flows. 

The presence of a cooling flow may act either to en- 
hance or reduce any metallicity gradient present in a cluster, 
with the details depending on how well the metals are mixed 
with the cooling gas. If the mixing is good then as the met- 
als are injected or deposited into the cluster environment 
they will primarily enrich the least dense X-ray gas at any 
radius, since this gas will have the largest volume filling fac- 
tor. This less-dense material will fiow furthest in towards the 
cluster center in the cooling fiow, whereas the denser, less- 
enriched material will be deposited at larger radii. Reiseneg- 
ger, Miralda-Escude & Waxman (1996) discuss such a model 
for the production of the metallicity gradient observed in 
the Centaurus Cluster, wherein the observed enrichment is 
primarily due to type la supernovae within the central clus- 
ter galaxy. A prediction of their model is that metallicity 
gradients should be most obvious in relatively low X-ray lu- 
minosity clusters, with higher ratios for their stellar/X-ray 
gas mass contents, which appears consistent with the detec- 
tions of strong gradients in the Centaurus and Virgo clusters, 
Abell 496 and AWM7. If the mixing of metals with the ICM 
is poor, however, localized regions of high metallicity ma- 
terial may cool faster than the surrounding gas, leading to 
the gradual reduction of any metallicity gradient apparent 
in the X-ray data, and to a decrease in the mean, emission- 
weighted metallicity of the cluster. Further discrimination 
on the origin of metallicity gradients in cooling fiow clusters 
will be possible using observations made with the Chandra 
Observatory and XMM-Newton. 



6 DEPROJECTION ANALYSIS OF THE 
ROSAT IMAGES 

The analysis of the ROSAT images has been carried out 
using an extensively updated version of the deprojection 
code of Fabian et al. (1981; see also White et al. 1997). 
Azimuthally-averaged X-ray surface brightness profiles were 
determined for each cluster from the ROSAT images. These 
profiles were background-subtracted, corrected for telescope 
vignetting and re-binned to provide sufficient counts in each 
annulus for the deprojection analysis to be extended beyond 
the cooling radii. 

With the X-ray surface brightness profiles as the pri- 
mary input, and under assumptions of spherical symmetry 
and hydrostatic equilibrium, the deprojection technique can 
be used to study the basic properties of the intracluster gas 
(temperature, density, pressure, cooling rate) as a function 
of radius. The deprojection code uses a monte-carlo method 
to determine the statistical uncertainties on the results and 
incorporates the appropriate HRI and PSPC spectral re- 
sponse matrices issued by GSFC. The cluster metallicities 
were fixed at the values determined from the spectral anal- 
ysis and the absorbing column densities were initially set to 
their respective Galactic values. 

The deprojection code requires the total mass profiles 
for the clusters, which define the pressure profiles, to be spec- 
ified. We have iteratively determined the mass profiles that 
result in deprojected temperature profiles (which approx- 
imate the mass-weighted temperature profiles in the clus- 
ters) that are isothermal within the regions probed by the 



ROSAT data and which are consistent with the best-fit tem- 
peratures determined from the ASCA spectra with spectral 
model E (or spectral models C and D for the Virgo and 
Coma clusters, respectively; Section 3). The validity of the 
assumption of isothermal mass-weighted temperature pro- 
files in the cluster cores is discussed by Allen (1998). The 
use of approximately isothermal deprojected temperature 
profiles also allows for a direct comparison between the mass 
deposition rates from the cooling flows determined from im- 
age deprojection and spectral analyses (see Section 7.1). 

The dominant components of the cluster mass profiles 
have been parameterized using the model for an isother- 
mal sphere given in Equation 4-125 of Binney & Tremaine 
1987) with adjustable core radii, re, and velocity disper- 
sions, o. The core radii were adjusted until the deviations 
from isothermality in the deprojected temperature profiles 
were minimized. The velocity dispersions were then adjusted 
until the temperatures determined from the deprojection 
code came into agreement with the values measured from 
the spectral analyses. An initial estimate for the pressure 
in the outermost radial bin used in the analysis is also re- 
quired by the deprojection code, although the uncertainties 
in these values do not significantly affect the results. 

For several of the clusters, these simple single- 
component mass models could not adequately satisfy the re- 
quirement for isothermal deprojected temperature profiles, 
particularly in the central regions of the clusters. However, 
for these systems it was also found that the introduction of 
a second 'linear' mass component (providing a constant con- 
tribution to the total mass per kpc radius c/ Thomas et al. 
1987) significantly improved the fits. These linear compo- 
nents (which may be intuitively associated with the central 
dominant galaxies in the clusters) were truncated at outer 
radii of 20 kpc. Their normalizations were adjusted itera- 
tively (together with the r^ and o values) in order to obtain 
the best match to the required isothermal temperature pro- 
files. 

The mass distributions determined from the deprojec- 
tion analysis are summarized in columns 3 — 5 of Table 9. 
The errors on the velocity dispersions are the range of val- 
ues that result in isothermal deprojected temperature pro- 
files that are consistent, at the 90 per cent confidence limit, 
with the spectrally-determined temperatures. These results 
may be used to provide direct comparisons with indepen- 
dent mass constraints from dynamical studies. We note that 
the mean core radius for the mass distributions in the cool- 
ing flow clusters of 48 ± 15 kpc is in good agreement with 
the value of ~ 50 kpc determined from the combined X-ray 
and gravitational lensing study of more distant, luminous 
lensing clusters by Allen (1998). 



7 THE PROPERTIES OF THE COOLING 
FLOWS 

In the following Section, we examine the basic properties 
of the cooling flows in the clusters. We compare the mass 
deposition rates determined from the spectral and image de- 
projection methods. We examine the evidence for intrinsic 
X-ray absorption in the clusters and compare our results 
with those from previous studies. Finally, we report detec- 
tions of extended infrared emission, spatially coincident with 
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the cooling flows, in several of the systems studied. The ob- 
served infrared fluxes and flux limits are shown to be in good 
agreement with the predicted values due to reprocessed X- 
ray emission from the cooling flows. 



7.1 The mass deposition rates 

The spectral and image deprojection analyses presented in 
Sections 3 and 6 provide essentially independent estimates 
of the mass deposition rates in the clusters. A comparison of 
the results obtained from these analyses therefore provides 
a flrm test of the validity of the cooling flow model. 

The deprojection method describes the X-ray emission 
from a cluster as arising from a series of concentric spherical 
shells. The luminosity in a particular shell, j, may be written 
as the sum of four components (Arnaud 1988). 



AMiH,+AMiA$i + 



E 



AMi{A<i>j + AHj 



(1) 



where AMj is the mass deposited in shell j, Hj is the en- 
thalpy of the gas in shell j, and A$j is the gravitational en- 
ergy released in crossing that shell. X/^-i ^-^^s is the mass 
flow rate through shell j, and AHj the change of enthalpy 
of the gas as it moves through that shell. The first term 
in equation nl thus accounts for the enthalpy of the gas de- 
posited in shell j. The second term is the gravitational work 
done on the gas deposited in shell j. The third and fourth 
terms respectively account for the gravitational work done 
on material fiowing through shell j to interior radii, and the 
enthalpy released by that material as it passes through the 
shell. 

In any particular shell, the densest material in the cool- 
ing flow is assumed to cool out and be deposited. Since the 
cooling time of this material will be short compared to the 
flow time, the cooling can be assumed to take place at a 
fixed radius. Thus, the luminosity contributed by the first 
term in equation 1 should have a spectrum appropriate for 
gas cooling at constant pressure from the ambient cluster 
temperature i. e. the same spectrum as the cooling fiow com- 
ponent in spectral model C (Section 3). For the bulk of the 
material continuing to flow inwards towards the cluster cen- 
tre, the cooling via X-ray emission is assumed to be offset 
by the gravitational work done on the gas as it moves in- 
wards. The emission accounted for in the second and third 
terms of equation 1 should therefore have a spectrum that 
can be approximated by an isothermal plasma at the ap- 
propriate ambient temperature for the cluster i.e. the spec- 
trum of the isothermal emission component also included in 
model C. Since the mass-weighted temperature profiles in 
the clusters are assumed to remain approximately isother- 
mal with radius, the luminosity contributed by the fourth 
term of equation nl should be relatively small. Overall, the 
integrated emission from the cooling flow can be expected to 
have a spectrum similar to that described by spectral model 
E. (In general, the results on the mass deposition rates deter- 
mined with spectral models C and E show good agreement; 
Table 4.) 

The mass deposition rates determined from the depro- 
jection analysis (Mi; Table 9) are the mass inflow rates mea- 
sured at the cooling radii in the clusters i.e. the radii at 



which the cooling time of the cluster gas first exceeds a Hub- 
ble time (1.3 x 10^" yr). Thus, if the cooling fiow model is 
correct, and if the cooling fiows have existed for a significant 
fraction of a Hubble time, the mass deposition rates deter- 
mined from the deprojection analysis should be compara- 
ble to the values measured independently from the spectral 
data. 

Fig. 4(a) shows the mass deposition rates determined 
from the spectral analysis (Ms) versus the results from the 
deprojection study (Mi). In all cases, except the Virgo Clus- 
ter, we have used the spectral results determined with model 
E as the ordinate, since this is the more physical of the two 
cooling fiow models and typically provides at least as good 
a fit to the data as the constant pressure approximation. 

The results plotted in Fig. 4(a) exhibit an approxi- 
mately linear correlation. A fit to the data with a power- 
law model of the form Ms ~ PM^, using the Akritas 
& Bershady (1996) bisector modification of the ordinary 
least-squares statistic (which accounts for the possibility 
of intrinsic scatter in the data) gives a best-fitting slope, 
Q — 1.15±0.33 and a normalization, P = 0.36±0.65 (where 
the errors are the standard deviations determined by boot- 
strap re-sampling). 

The spectral analysis presented in Section 3 also indi- 
cates that the cooling gas is intrinsically absorbed by equiva- 
lent hydrogen column densities of, typically, a few 10'^'^ atom 
cm~^. In principle, therefore, the mass deposition rates de- 
termined from the deprojection analysis should also be cor- 
rected for the effects of this absorbing material. Such cor- 
rections, however, depend upon both the spatial distribution 
and nature of the absorbing material, which at present re- 
main uncertain. Nevertheless, under the specific assumption 
that the intrinsic absorption is due to cold gas which lies 
in a uniform screen in front of the cooling flows, the effects 
of intrinsic absorption on the deprojection results may be 
estimated. 

The corrections for intrinsic absorption due to a uni- 
form foreground screen of cold gas have been carried out 
by re-running the deprojection code with the absorbing col- 
umn densities set to the total values determined with spec- 
tral model E (model C for the Virgo Cluster). Since the 
intrinsic column densities measured with model E are red- 
shifted quantities, we have set the total column densities 
used in the revised deprojection analysis, which assumes 
zero redshift for the absorber, to values A''h-I~ AA'^h/(1 + -z)"^, 
where Nn is the Galactic column density.) The absorption- 
corrected mass deposition rates (Mc) so determined are 
summarized in Table 10. Fig. 4(b) compares the results on 
the mass deposition rates from the spectral analysis with the 
absorption-corrected deprojection values. A fit to the data 
with a power-law model of the form Ms = PM^j gives a 
slope, g = 1.09 ±0.32 and a normalization, P = 0.18 ±0.37 
(where the errors are the standard deviations determined by 
bootstrap re-sampling). 

We conclude that the results on the mass deposition 
rates determined from the spectral and imaging methods 
exhibit an approximately linear correlation and marginal 
agreement in normalization in the case where no corrections 
for intrinsic absorption on the imaging results are made. 
The application of simple correction factors, appropriate for 
absorption by a uniform screen of cold gas, lead to the imag- 
ing results exceeding the spectral values by typical factors 
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of ~ 2 — 4, although the details of these corrections remain 
highly uncertain at present. As we shall discuss in Section 
8, an improved agreement between the imaging and spectral 
results on the mass deposition rates is obtained once the dif- 
ferent ages of the cooling flows are also accounted for in the 
analysis. 



and that the absorbing material lies at zero redshift. Model 
C has been used in a number of previous studies {e.g. White 
et al. 1991; Allen & Fabian 1997) and is included here for 
comparison purposes. The results obtained for the present 
sample of clusters using model C are also listed in Table 11. 
The mean excess column density for the cooling-flow clusters 



7.2 Intrinsic X-ray absorption 

The negative residuals at low energies [E < O.SkcV) in Fig. 
2 provide a clear indication of the presence of excess absorp- 
tion, over and above the nominal Galactic values (Dickey & 
Lockman 1990) along the lines of sight to the clusters. The 
analyses of the full ASCA data sets with spectral models 
B, C, D and E also indicate the presence of intrinsic X-ray 
absorbing material in the clusters. However, the precise col- 
umn density measurements, summarized in Table 11, are 
sensitive to the spectral model used. 

Using the simple isothermal model with free-fitting ab- 
sorption (model B), we measure a mean excess column den- 
sity for the seven cooling-flow clusters of 2.4 ± 2.8 x 10^" 
atom cm~^. We note that the value measured for the Perseus 
Cluster (Abell 426) is negative, supporting the conclusion 
drawn in Section 3.3 that the nominal Galactic column den- 
sity determined from HI studies (A'^h ~ 1.49 x 10 atom 
cm~^; Dickey & Lockman 1990) overestimates the true col- 
umn density along the line of sight to the cluster core. Inter- 
estingly, the excess column density measured for the Coma 
cluster (1.5 ± 0.5 x lO^'' atom cm~^) is consistent with the 
mean value determined for the cooling-flow clusters, using 
the same spectral model. Our results on the excess column 
densities measured with spectral model B may be compared 
to the 90 per cent confldence upper limits on the size of 
any systematic excess column density present in the SO and 
SI detectors due to calibration errors of < 8 x 10^^ and 
< 2 X 10^" atom cm^^, respectively (Orr et al. 1988; these 
values are valid for observations made before mid 1996. See 
also Dotani et al. 1996). 

The results on the intrinsic column densities for the 
cooling-flow clusters are quite different, however, when the 
more-sophisticated multiphase spectral models (C, D and E) 
are used. As discussed in Section 3, models C, D and E gen- 
erally provide a significantly better fit to the ASCA data for 
the cooling-flow clusters than the simple single-phase mod- 
els (A and B). Using spectral models C and E, which are our 
preferred models in that they provide a consistent descrip- 
tion for the spectral and imaging X-ray data for the clusters 
(Sections 7.1, 8.3), we measure a mean intrinsic column den- 
sity acting on the cooling-flow components of 3 — 4 x lO'^^ 
atom cm~^. Using the two-temperature model (model D) we 
determine a similar mean intrinsic column density. Thus, the 
typical intrinsic column densities measured with the multi- 
phase spectral models are approximately an order of mag- 
nitude larger than those determined using the single-phase 
model B. These results are in good agreement with those 
presented by Allen (2000) from an analysis of larger sample 
of more luminous, distant cooling flow clusters. 

We have also measured the intrinsic column densities 
in the clusters using one further spectral model, which we 
refer to as model C. Model C is identical to model C except 
that it assumes that the intrinsic absorption acts on the 
entire cluster spectrum, rather than just the cooling gas. 



Allen & Fabian (1997) presented results from an X-ray 
colour deprojection study of 18 clusters observed with the 
ROSAT PSPC, from which they determined intrinsic col- 
umn densities across the central 30 arcsec (radius) regions 
of Abell 426, 496, 1795, 2199 and 2597 (using spectral model 
C) of 0.50 ± 0.55, 2.57 ± 0.62, 2.20 ± 0.20, 2.67 ± 0.25 and 
2.36 ± 0.58 X lO^'' atom cm^^, respectively. The results for 
Abell 496, 1795 and 2199 are in reasonable agreement with 
those presented here, although our value for Abell 2597 is 
~ 3 times larger than the Allen & Fabian (1997) result. We 
note, however, that the ASCA observations of Abell 2597 
were carried out in 1996 December and are likely to have 
been affected by radiation damage to the SIS detectors, 
which is known to lead to overestimates of the measured 
column densities in observations made after (approximately) 
late 1995 (Hwang et al. 1999). The Allen & Fabian (1997) 
result for Abell 426 was determined using a lower value for 
the Galactic column density of 1.37 x 10^^ atom cm^^ (Stark 
et al. 1992). If we also adopt this lower value, our ASCA 
measurement for the intrinsic column density in this clus- 
ter, using model C, rises to O.IOIq^q x 10^°, consistent with 
the previous PSPC result. 

Six of the cooling-flow clusters included in the present 
work were also studied by White et al. (1991), using Ein- 
stein Observatory SSS data. The analysis of the SSS obser- 
vations (which covered central, three arcmin radius circular 
apertures in the clusters) was also carried out using spectral 
model C. White et al. (1991) presented measurements of 
the intrinsic column densities in AbeU 426, 496, 1795, 2199 
and the Virgo cluster of 13 ± 2, 20+*, 8 ± 3, 14 ± 2 and 
15±2 X 10^° atom cm~^, respectively. (The SSS observation 
of the Centaurus cluster was offset from the centroid listed 
in Table 3 by ~ 4 arcmin, prohibiting any direct compari- 
son with the present work.) The intrinsic column densities 
measured by White et al. (1991) are significantly larger (by 
factors ranging from two to > an order of magnitude) than 
the values determined from the ASCA and ROSAT data. 
The ASCA result for the non-cooling fiow Coma cluster, 
measured with spectral model B, of AA''h ~ 1.5 ± 0.5 x lO'^" 
atom cm^^ is consistent with the White et al. (1991) limit 
of AiVn < 5 X 10^° atom cm'^ 

Finally, we have examined the constraints that may be 
placed on the covering fraction, /, of the intrinsic X-ray ab- 
sorbing material using the ASCA data. Using spectral model 
E, we measure a typical best-fit covering fraction of unity 
and can rule out covering fractions significantly less than 90 
per cent. The exception to this is the result for the Virgo 
cluster, for which we measure a best-fit covering fraction 
with spectral model E, / = 0.66+q ;[g (for spectral model C 
we find / = 0.86+QgY). The results on the covering fractions, 
and the agreement of the ASCA and PSPC results, suggest 
that the White et al. (1991) results are likely to have sys- 
tematically over-estimated the intrinsic column densities in 
clusters, although the qualitative result on the detection of 
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significant intrinsic X-ray absorption in cooling-flow clusters 
is confirmed. 



7.3 The luminosity reprocessed in other 
wavebands 

Tlie luminosities absorbed at X-ray wavelengths must even- 
tually be reprocessed in other wavebands. If the absorbing 
material is dusty (as is likely to be the case in the central 
regions of the clusters e.g. Voit & Donahue 1995; Fabian 
et al. 1994a) then the bulk of this reprocessed emission is 
likely to emerge in the far infrared. Table 12 summarizes the 
reprocessed luminosities (i.e. the bolometric luminosities ab- 
sorbed within the clusters) measured with spectral model E 
(or model C for the Virgo Cluster). The reprocessed lumi- 
nosities range from ~ 4 x 10*^ ergs^^ for the Virgo Cluster 
to ~ 5 X lO""* ergs"^ for AbeU 2597. 

Table 12 also lists the observed 60 and lOO/im fluxes 
within a four arcmin (radius) aperture centred on the X-ray 
centroids for the clusters (Table 3) . The observed fluxes were 
measured using the IPAC SCANPI software and co-added 
IRAS scans. (The median of the co-added scans was used.) 
The root-mean-square deviations in the residuals external to 
the source extraction regions after baseline subtraction, and 
the in-scan separations (in arcmin) between the peaks of the 
100/im emission and the X-ray centres, are also included in 
the table. 

Three of the clusters hsted in Table 12, AbeU 426, the 
Centaurus Cluster and AbeU 2199, provide convincing detec- 
tions of IQQixra flux associated with the peaks of their X-ray 
emission. For the Centaurus cluster and AbeU 2199 (and the 
slightly offset source associated with AbeU 496), the infrared 
emission appears spatially extended. For AbeU 426, the de- 
tected infrared flux is point-like and presumably associated 
with the active nucleus of the central galaxy, NGC1275. Cox, 
Bregman & Schombert (1995) have previously reported de- 
tections of infrared emission from the central regions of the 
Centaurus Cluster and AbeU 2199, from an analysis of the 
same IRAS scan data. These authors did not detect sig- 
nificant (> 98 per cent confidence) infrared emission from 
within 1 arcmin (radius) of the dominant galaxies in AbeU 
496, 1795 or 2597, which is consistent with our results. The 
presence of a strong, point-like infrared source in AbeU 426 
was previously reported by Wise et al. (1993) who measured 
60 and lOO/im fluxes associated with this source of 6.48 and 
8.73 Jy, respectively. 

Following Helou et al. (1988) and Wise et al. (1993), the 
total infrared luminosities in the 1 — lOOO/im band may be 
estimated from the observed IRAS fluxes using the relation 

Li-ioooMm ~ 2.8 X 10'*''(^)2(2.58S'6o + Sioo)crgs-\ (2) 

where Sso and 6*100 are the 60 and lOO/im IRAS fluxes in 
units of Jy. This relation assumes a dust temperature of 
~ 30K and an emissivity index n in the range — 2, where 
the emissivity is proportional to the frequency, v" . We as- 
sociate a systematic uncertainty of ±30 per cent with the 
estimated 1 — lOOO/im luminosities, which is combined in 
quadrature with the random errors. (In most cases, the sys- 
tematic error exceeds the random errors. Note that the sys- 
tematic errors associated with the measured IRAS fluxes 
may be ^ 30 per cent in some cases e.g. see the note on 



the Virgo Cluster below.) The 1 — 1000/u,m luminosities cal- 
culated from this relation are summarized in Table 12. Fig. 
5 shows the results plotted as a function of the predicted 
reprocessed luminosities. 

For the Centaurus Cluster and AbeU 2199, the esti- 
mated 1 — 1000/im luminosities are in good agreement with 
the reprocessed luminosities predicted from the X-ray mod- 
els. The observed Seo/^ioo flux ratios for these clusters in- 
dicate dust temperatures in the range 20 — 50K. For AbeU 
426, the 1 — 1000/im luminosity exceeds the predicted repro- 
cessed X-ray luminosity by a factor ~ 4, consistent with the 
presence of a strong, intrinsically absorbed active nucleus 
in this system (c/ Section 4). For AbeU 496, the observed 
flux also exceeds the predicted flux by a factor of ~ 4, al- 
though the ~ 3 arcmin spatial offset between the 100/im 
and X-ray centroids suggests that the infrared emission is 
likely to originate, at least in part, from some source other 
than the cooling flow. (This is also evidenced by the im- 
plausibly low 56o/5ioo flux ratio measured for the cluster). 
For the Virgo Cluster, the predicted reprocessed X-ray lumi- 
nosity exceeds the 1 — 1000/im luminosity determined from 
the IRAS data by a factor of three. However, this is un- 
surprising since the SCANPI software is not well suited to 
measuring fluxes from highly extended objects (the cooling 
flow in the Virgo Cluster has the largest angular extent of 
any of the systems studied in this paper, although this may 
also affect the IRAS results for the Perseus Cluster.) Finally, 
we note that it if the bulk of the absorbing material in the 
clusters is very cold (T < lOK), as suggested by some mod- 
els {e.g. Johnstone, Fabian & Taylor 1998 and references 
therein), then a signiflcant fraction of reprocessed emission 
may emerge in the sub-mm band. 



8 MEASURING THE AGES OF COOLING 
FLOWS 

The natural state for a regular, relaxed cluster of galaxies 
appears to include the presence of a cooling flow in its core 
{e.g. Edge et al. 1992; Peres et al. 1998). Simulations sug- 
gest that cooling flows are only likely to be disrupted to the 
extent that they are 'turned off' when a cluster merges with 
a subcluster of comparable size {e.g. McGlynn & Fabian 
1984). In this section, we discuss three methods which may 
be used to measure the ages of cooling flows from the X-ray 
data and apply them to the ASCA and ROSAT observations 
of the nearby cluster cooling flows. 



8.1 Method 1: X-ray colour deprojection 

The flrst constraints on the ages of cooling flows in clus- 
ters were presented by Allen & Fabian (1997) from their 
X-ray colour deprojection study of 18 nearby systems ob- 
served with the ROSAT PSPC. These authors used their 
X-ray colour deprojection technique to measure the spa- 
tial distributions of cool(ing) gas in the clusters and com- 
pared their results with the predictions from simple, time- 
dependent cooling flow models. Essentially, Allen & Fabian 
(1997) identified the ages of the cooling flows with the cool- 
ing time of the ICM at the outermost radii at which signifi- 
cant cooling is observed. These authors provide cooling-fiow 
age measurements for seven clusters, four of which are also 
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included in the present work. Their results for AbcU 496, 
1795, 2199 and 2597 are 3.3 - 6.8, 3.6 - 6.7, 5.3 - 7.7 and 
2.5 — 5.1 Gyr, respectively. 

8.2 Method 2: Comparison of the spectral and 
imaging mass deposition rates 

A second estimate of the ages of cooling flows may be ob- 
tained by comparing the essentially independent results on 
the mass deposition rates determined from the spectral and 
image deprojection studies. The deprojection analysis, dis- 
cussed in Section 6, measures the mass deposition profiles 
in the clusters i. e. the mass deposition rates as a function of 
radius. The integrated mass deposition rates inferred from 
the deprojection analysis are the mass deposition rates inte- 
grated out to some critical radius (rcooi) at which the mean 
cooling time of the X-ray gas (icooi) becomes equal to an as- 
sumed age for the system. (The results presented in Table 9 
are for an assumed age of 1.3 x 10^° yr). The spectral data, in 
contrast, provide a measure of the current, total mass depo- 
sition rates in the clusters. Thus, the ages of the cooling flows 
may simply be associated with the cooling time of the X-ray 
gas at the radii where the integrated mass deposition rates 
determined from the deprojection analysis become equal to 
the values measured directly from the spectral data. 

The primary uncertainty with this method lies in the 
corrections for the effects on intrinsic absorption, especially 
on the deprojection results. The larger the correction for in- 
trinsic absorption at a particular radius, the larger the mass 
deposition rate inferred at that radius and, therefore, the 
smaller the age of the cooling flow. The cooling-flow ages 
determined with this method, both with and without cor- 
rections for intrinsic absorption due to a uniform foreground 
screen of cold gas on the deprojection analysis, are summa- 
rized in columns 3 and 4 of Table 13. 

8.3 Method 3: Breaks in the deprojected mass 
deposition profiles 

The X-ray colour deprojection method (method 1) identifles 
the age of a cooling flow with the mean cooling time of the 
cluster gas at the outermost radius where the spectral signa- 
ture of cooling gas is observed. Since the presence of cooling 
gas will tend to enhance the X-ray luminosity of a cluster, 
the outermost radius at which cooling occurs may also be 
expected to be associated with a 'break' in the X-ray surface 
brightness profile and, more evidently, the mass deposition 
profile determined from the deprojection analysis. 

We have searched for the presence of breaks in the mass 
deposition profiles determined from the deprojection analy- 
sis by fitting these profiles with broken power-law models. 
Only those data from radii interior to the 90 percentile upper 
limits to the cooling radii were included in the modelling, 
with the exception of the Virgo cluster where, due to system- 
atic uncertainties associated with the effects of the detector 
support structure, the analysis was limited to the central 
78 kpc (15 arcmin) radius. An unweighted, ordinary least- 
squares statistic was used in the fits. (A reduced x^ value of 
1.0 was assumed in estimating the Ax^ = 2.71 errors.) The 
results from the fits are summarized in Table 14. 

Several of the clusters studied {e.g. Abell 426, 496) ex- 
hibit clear breaks in their mass deposition profiles at radii 



r < Tcooi. In other systems, the identification of a break 
radius is less clear, although still formally significant (e.g. 
AbcU 2199). Following method 1, we identify the ages of the 
cooling fiows with the mean cooling time of the cluster gas at 
the break radii. In calculating the ages, we assume that the 
X-ray emission at the break radii is unaffected by intrinsic 
absorption, which is reasonable if the absorbing material in 
the clusters is gradually accumulated by the cooling flows. 

The only cluster for which we do not detect a significant 
break in the mass deposition profile is the Centaurus Clus- 
ter which, interestingly, also has the largest age measured 
with method 2. (The age measurement for the Centaurus 
Cluster with method 2 is consistent with the cooling flow 
having survived intact since the cluster was formed. In this 
case no break in the mass deposition profile is expected). 
Fig. 6 shows the mass deposition profiles for a representa- 
tive subsample of the clusters studied, with the best-fitting 
broken power-law models overlaid. The mean cooling time 
of the cluster gas as a function of radius is shown in the 
lower panels. 

We have measured the slopes of the mass deposition 
profiles interior and external to the break radii in the clus- 
ters. (These slopes may be used to constrain the range of 
plausible isothermal cooling flow models discussed in Sec- 
tion 3.2). The mean slope interior to the break radii (i.e. in 
the regions where mass deposition is expected to occur) is 
1.37 ± 0.27. Beyond this point, the mean slope flattens to a 
value of 0.59 ± 0.21. 

Finally, in Table 14, we also list the integrated mass 
deposition rates interior to the break radii, both corrected 
and uncorrected for absorption due to a uniform foreground 
screen of cold gas (as described in Section 7.1). Fig. 7 shows 
these values plotted as a function of the mass deposition 
rates determined from the spectral data. (The Centaurus 
Cluster does not exhibit a clear break within the cooling ra- 
dius and we assume Merk = Mc or Merk = Mi from Table 
10, as appropriate.) In general, the agreement between the 
results is good, particularly for the uncorrected data. (Re- 
call that significant uncertainties remain in the corrections 
for intrinsic absorption, due both to the unknown geome- 
try and physical nature of the absorbing matter. See also 
Section 9.1). A fit to the results with a power-law model 
of the form Ms — PM^^^^, using the Akritas & Bershady 
(1996) bisector modification of the ordinary least-squares 
statistic, gives a slope Q = 1.08 ± 0.14 and a normalization, 
P — 0.33 ± 0.25 for the absorption-corrected data, and a 
slope Q = 1.11 ± 0.19 and a normalization, P — 0.70 ± 0.61 
for the uncorrected results (where the errors are the stan- 
dard deviations determined by bootstrap re-sampling). 



8.4 Comparison of results 

Table 13 summarizes the ages for the cooling flows deter- 
mined using the three different methods. These results are 
also shown in graphical form in Fig. 8, where we have av- 
eraged the absorption-corrected and uncorrected results de- 
termined with method 2. In general, the results on the ages 
show good agreement, suggesting that the individual meth- 
ods may indeed be used to reliably measure the ages of the 
cooling flows. 

One source of uncertainty with the age measurements is 
the relation between the cooling time {i.e. the time taken for 
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the gas to cool to zero degrees at constant pressure) of the 
ICM at the outermost radii where cooUng gas is observed 
and the ages of the cooling flows, fn detail, the outermost 
radius at which significant mass deposition occurs will de- 
pend upon the spectrum of density inhomogeneities in the 
cluster gas, which at present remains unknown. 

As discussed above, the effects of intrinsic absorption 
provide an additional source of uncertainty. In particular, 
the ages determined with method 2 (the comparison of the 
spectral and deprojected mass deposition rates) are sensi- 
tive to such effects, although the results based on the break 
radii (method 3) are insensitive to any uniform change in 
column density across the cooling flows, and those from the 
X-ray colour deprojection analysis (method 1) are largely 
independent of column density uncertainties (other than in 
the determination of the cooling time at the outermost edge 
of the cooling flows, where we assume zero intrinsic absorp- 
tion). Future observations with the Chandra Observatory 
and XMM-Newton will directly address such uncertainties, 
permitting spatially-resolved spectroscopic determinations 
of both the mass deposition and intrinsic absorption profiles 
in clusters on spatial scales > a few arcsec. 

Taking into account the uncertainties associated with 
the different methods, we identify method 1 as the most ro- 
bust, followed by method 3, with method 2 being the most 
sensitive to systematic effects. Method 1 requires the best 
data, with good spatial and moderate spectral resolution 
and a high signal-to-noise ratio. (To date this method has 
has only been applicable to PSPC observations of a few 
of the nearest, brightest clusters of galaxies). In contrast, 
method 3 should be applicable to a large number of ROSAT 
HRI observations. We note that method 2 is, at present, 
also complicated by the fact that it combines results from 
different detectors, which introduces additional systematic 
uncertainties. 

Averaging the absorption-corrected and uncorrected re- 
sults determined with method 2 (which provides a reason- 
able estimate of our uncertainty in these quantities) and 
taking an unweighted mean of the results determined with 
the three separate methods, we arrive at the following ages 
for the cooling flows in the clusters: Abell 426 (2.5 Gyr), 
AbeU 496 (5.4 Gyr), the Virgo Cluster (2.6 Gyr), the Cen- 
taurus Cluster (7.1 Gyr), AbeU 1795 (5.4 Gyr), AbeU 2199 
(5.9 Gyr) and AbeU 2597 (3.2 Gyr). 



9 DISCUSSION 

9.1 The possible role of heating processes 

The age-corrected mass deposition rates determined from 
the image deprojection analysis, with no corrections for the 
effects of intrinsic absorption applied, are generally in good 
agreement with the spectrally-determined values. However, 
the absorption-corrected, age-corrected deprojection results 
are typically about a factor of two larger than the spec- 
tral measurements. Although this discrepancy may simply 
be due to incorrect assumptions about the spatial distribu- 
tion and physical nature of the absorbing medium (Voit & 
Donahue 1995, Arnaud & Mushotzky 1998 and Allen 2000 
discuss the possibility that the intrinsic absorption in clus- 
ters may be primarily due to dust grains) it is sensible to 
also consider other possibilities. 



The location of the break radii in the clusters is cru- 
cial in measuring the integrated mass deposition rates from 
the imaging data. One possibility for reducing the imaging 
values is to associate the ages of the cooling flows with the 
cooling time of the cluster gas at radii of about one half of 
the observed break radii. This would, however, lead to sig- 
nificantly shorter ages for the cooling flows, which would be 
puzzling for a phenomenon which appears to be so common. 
Moreover the X-ray colour deprojection technique (Method 
1) clearly shows that cool(ing) gas extends out to about the 
break radii. The assumption of constant metallicity within 
the cooling flows (where significant metallicity gradients are 
observed; Section 5.3) will also introduce systematic errors, 
although these are unlikely to result in discrepancies as large 
as a factor of two. The use of a larger value for rj in the 
isothermal cooling flow spectral models (Table 5) can also 
slightly reduce the discrepancy between the spectral and 
absorption-corrected imaging mass deposition rates, but not 
fully account for it. It remains possible that a combination 
of some or all of these factors, together with issues such 
as the correctness of the plasma code, details depending on 
the assumed element abundances (the current cooling flow 
models do not account for variable element abundance ra- 
tios) and/or the geometry and effects of the magnetic fleld, 
are responsible. 

A further plausible alternative is that there may be 
some heating as well as cooling in the central region. Many 
authors have proposed various heating models (see Fabian 
1994 for references) but none has been found to be com- 
pelling and essential in the face of the clear spectral evidence 
for cool, and probably cooling, gas components. Nevertheless 
the presence of radio sources at the centres of most cooling 
flows demonstrates that at least some heating occurs. An 
important point with any viable heating process, given the 
X-ray observations, is that it must allow some of the gas to 
continue cooling, which means that the heating rate must 
be say one third to one half of the total cooling rate. Such 
a close balance suggests feedback i.e. the heat source is fed 
by the cooling gas. 

One obvious solution is that the central engines in the 
radio sources are fuelled by the flows (Bailey 1980; Nulsen, 
Stewart & Fabian 1984; Fabian & Crawford 1990). This has 
the useful properties that the radio sources are spatially ex- 
tended and can therefore supply heat (probably via low en- 
ergy cosmic rays, depending on the magnetic topology and 
connection to the various gas phases, or by shocks) over an 
extended region, and that the cooling gas fueling the source 
will, in a multiphase flow, be the hottest phase in which the 
effects of cosmic ray and/or shock heating are greatest (in 
general, the hottest phase will have the longest cooling time 
and shock or cosmic ray heating will increase its cooling 
time by the greatest absolute amount. Presumably heating 
will also be strongest near the centre). This hottest phase 
is then likely to be the phase through which the feedback 
operates. If there is temporarily too much heating, this may 
stop the cooling and reduce the central fuel supply to the 
radio source. Conversely, if there is too little heating the ac- 
cretion rate onto the nucleus will rise (for a discussion in the 
context of single-phase gas models see e.g. Tucker & David 
1997). 

Pedlar et al. al (1992) and Owen & Eilek (1998) have 
emphasised the importance of heating by the central radio 



© 0000 RAS, MNRAS 000, 000-000 



ASCA and ROSAT observations of nearby cluster cooling flows 15 



sources in the Perseus Cluster and Abell 2199, respectively. 
Our results provide an indication of the possible strength of 
this heating in such systems. The heating may occur inter- 
mittently, but overall it appears that it may cause up to a 
factor of two reduction in the total mass deposition rates 
deduced from the image deprojection method, once absorp- 
tion corrections have been made. We note that these rates 
remain close to the values obtained when no absorption cor- 
rection is made and are therefore similar to the pre- 1990 
published values, which have always been regarded as un- 
certain by about a factor of two {e.g. Arnaud 1988; Fabian 
1994) due to uncertainties in the gravitational potentials of 
the clusters. 

The tentative identification of a heating signature dis- 
cussed here does not affect the question of what becomes of 
the cooled material deposited by cooling flows. The spectral 
evidence from ASCA and other X-ray data suggest that tens 
to hundreds of solar masses of gas cool out of cooling flows 
per year. Over the billion year age of a flow this amounts 
to 5 X 10^" to greater than 10^^ M0 of cooled material, of 
which only a small fraction is directly observed at other 
wavelengths. 



9.2 The masses of merging subclusters that 
disrupt cooling flows 

Observations and theoretical models suggest that cooling 
flows are likely to be disrupted when a cluster merges with 
a subcluster of comparable size {e.g. McClynn & Fabian 
1984; Edge, Stewart & Fabian 1992; Allen 1998). When the 
density in the core of the merging subcluster is sufficiently 
high, the gas in the central regions of the dominant cluster 
is likely to be shocked and disrupted as the systems merge, 
leading to complex X-ray morphologies, temperature inho- 
mogeneities and an inflation of the X-ray core radii (Navarro, 
Frenk & White 1995; Schindler 1996; Roettiger, Burns & Lo- 
ken 1996). The time taken for the cooling flow to reform after 
the dark matter distributions merge and hydrostatic equilib- 
rium is restored will depend upon the degree to which cool 
gas originally present in the central regions of the clusters 
is heated and mixed with the surrounding ICM. However, 
the fact that ~ 70 per cent of bright, nearby clusters con- 
tain cooling flows (Peres et al. 1998) suggests that cooling 
flows are 'turned off' on average for ~ 30 per cent of the 
time. In this case, we may expect the typical time taken for 
a cooling flow to reestablish itself following a major subclus- 
ter merger event to be of the order of a few Gyr. Given the 
results on the ages of the cooling flows presented in Section 
8, we may therefore expect the present sample of clusters to 
have experienced their last major subcluster merger events 
during the last 5 — 10 Gyr (with the Coma, Perseus and 
Virgo clusters having experienced the most recent merger 
events and the cooling flow in the Centaurus Cluster having 
been undisturbed for the longest period). 

Assuming that subcluster mergers are the primary 
cause of the disruption of cooling flows, we can estimate 
the masses of the merging subclusters which cause signif- 
icant disruption by integrating the merger rates predicted 
by hierarchical formation models, using parameters appro- 
priate for rich clusters. Such rates are given by the extended 
Press-Schechter formalism of Lacey & Cole (1993). Exam- 



ples of their application to rich clusters are given by Cava- 
liere, Menci & Tozzi (1999). 

The results for a cluster mass of 10^^ Mq are shown 
in Fig. 9, where we plot the fraction of clusters that have 
undergone a merger with a subcluster of mass greater than 
10 — 80 per cent of the final cluster mass, within a given 
interval. (We assume O = 1.0, A = 0, 6 = 1 and a cosmic 
fluctuation index n — —1.5.) It is clear from the figure that 
for there to have been a good chance (> 50 per cent) of a 
merger within the last 10 Gyr, subclusters with masses of 
< 30 — 40 per cent of the total final mass are required. Major 
mergers involving subclusters of similar mass {i.e. half of the 
final mass) have a probability of ~ 20 per cent within that 
time. The results are similar for an open cosmology with 
Q — 0.3 (other cosmological parameters the same) for which 
subclusters with masses of < 20 — 30 per cent of the total, 
final mass have an > 50 per cent chance of merging within 
the last 10 Gyr and the probability of a merger with a sub- 
cluster of similar mass is ~ 10 per cent. Thus, if subcluster 
mergers are the primary cause of the disruption of cooling 
fiows in rich clusters, it seems likely that subclusters with 
masses of ~ 30 per cent of the final cluster mass can disrupt 
cooling fiows. 

The results for a cluster of mass lO^** M0 (a mass sim- 
ilar to that of the lowest luminosity clusters in our sample) 
show a much higher rate of significant mergers, with a ~ 50 
per cent chance of a major merger within the last 5 Gyr. 
It is therefore not clear why the ages of the cooling fiows 
in the present sample of objects appear roughly indepen- 
dent of their luminosity (and therefore mass). A possible 
explanation is that mergers may not be the only cause of 
cooling flow disruption. Activity in the nuclei of the central, 
dominant galaxies may also be important. If the ages of the 
cooling flows (taking into account the time required for the 
flows to settle down) are around 8 Gyr, then the disruption 
is likely to have occurred at redshift of 2; = 1 — 2 and so 
could be related to the peak of quasar/radio activity. 

The application of similar techniques to future obser- 
vations made with the Chandra Observatory and XMM- 
Newton should permit precise measurements of the ages of 
cooling flows for large samples of clusters within redshifts 
z<l. 



10 CONCLUSIONS 

The main conclusions that may be drawn from this paper 
may be summarized as follows: 

(i) We have demonstrated the need for multiphase mod- 
els to consistently explain the spectral and imaging X-ray 
data for the nearest, brightest cooling-flow clusters. In gen- 
eral, the diffuse, thermal emission from these systems can 
be modelled using either two temperature or cooling flow 
models, with excess absorption acting on the cooler emis- 
sion components. 

(ii) The mass deposition rates for the cooling flows, in- 
dependently inferred from the ASCA spectra and ROSAT 
images, exhibit reasonable agreement, providing support for 
the basic cooling flow picture. We have discussed the sys- 
tematic uncertainties associated with the results, which are 
present at the factor ~ 2 level, and commented on the pos- 
sible effects of intermittent heating in the cooling flows. 
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(iii) The ASCA spectra for the central regions of the 
Virgo and Centaurus clusters place firm constraints on the 
abundances of iron, magnesium, silicon and sulphur in the 
ICM. Our results firmly favour enrichment models in which 
a high mass fraction (70 — 80 per cent) of the iron in the 
ICM in the central regions of these clusters (in which strong 
metallicity gradients are observed) is due to Type la super- 
novae. 

(iv) We have confirmed the presence of intrinsic X-ray 
absorption in the ASCA spectra for the clusters using a 
variety of spectral models. The measured column densities 
are sensitive to the spectral models used and range (on av- 
erage) from a few 10^° atom cm~^ using a simple single- 
temperature model to a few lO^'^ atom cm~^ using our pre- 
ferred cooling-fiow models. 

(v) We have reported detections of 60 and lOO/im in- 
frared emission, coincident with the X-ray centroids, from a 
number of the systems studied (confirming the original de- 
tections by Wise et al. 1993 and Cox et al. 1995). For the 
Centaurus Cluster and Abell 2199, the coadded IRAS scan 
data show that the infrared emission is spatially extended. 
The infrared fiuxes and fiux limits are in good agreement 
with the predicted values due to reprocessed X-ray emission 
from the intrinsic X-ray absorbing material in the cooling 
flows. 

(vi) We have discussed three methods which may be 
used to measure the ages of cooling flows from the X-ray 
data. The application of these techniques to the present sam- 
ple of objects indicates cooling-fiow ages of between 2.5 and 
7 Gyr. If the ages of cooling fiows in the most massive clus- 
ters clusters are primarily set by subcluster merger events, 
then our results suggest that mergers with subclusters with 
masses of ~ 30 per cent of the final cluster mass are sufficient 
to disrupt cooling fiows. 
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Table 1. Summary of the ASCA Observations. Column 2 lists the optically-determined redshifts for the clusters (for the Virgo Cluster 
we assume a distance of 18 Mpc). Column 3 summarizes the dates of the observations. Columns 4 — 7 list the net exposure times (in 
seconds) in each of the four detectors after all screening and cleaning procedures were carried out. 



Cluster 


z 


Date 


SO 


SI 


G2 


G3 


Abell 426 


0.0183 


1993 Aug 06 


13910 


6814 


15618 


14916 


Abell 496 


0.0320 


1993 Sep 21 


29685 


21515 


40050 


39986 


Virgo 


18Mpc 


1993 Jun 07 


13633 


14668 


16874 


16874 


Centaurus 


0.0104 


1995 May 19 


68820 


67650 


70981 


70999 


Coma 


0.0232 


1993 Jun 14 


8732 


6637 


9312 


9312 


Abell 1795 


0.0634 


1993 Jun 16 


31284 




37649 


37641 


Abell 2199 


0.0309 


1993 Jul 25 


19153 


13932 


29758 


29758 


Abell 2597 


0.0852 


1996 Dec 06 


33933 


34697 


39249 


39239 



Table 2. The radii of the circular extraction regions used in the analysis of the ASCA SIS data (in arcmin and kpc) and the chip modes 
used in the observations (either 1,2 or 4-CCD mode). The numbers in parentheses indicate the number of chips contributing to the 
extracted spectra. For the CIS data a fixed extraction radius of 6 arcmin was used. 



Cluster 


SO 


SI 


SIS Chip Mode 




(amin/kpc) 


(amin/kpc) 




Abell 426 


4.9/152 


4.7/145 


4(1) 


Abell 496 


3.9/206 


3.2/169 


4(1) 


Virgo 


4.9/25.5 


4.1/21.3 


4(2) 


Centaurus 


— 


3.7/66.0 


1(1) 


Coma 


10.0/389 


9.5/370 


4(4) 


Abell 1795 


4.5/447 


— 


4(2) 


Abell 2199 


5.3/271 


4.2/215 


4(2) 


Abell 2597 


4.5/580 


3.7/477 


2(2) 
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Table 3. Summary of the ROSAT Observations. Columns 2 and 3 list the dates of the observations and the instrument used. Column 4 
lists the exposure times (in seconds). Where more than a single observation of a cluster was made, details for each observation are given. 
Columns 5 and 6 list the coordinates for the centroids of the X-ray emission from the clusters. 



Cluster 


Date 


Instrument 


Exposure 


R.A. (J2000). 


Dec. (J2000.) 


Abell 426 


1992 Feb 02 


PSPC 


4787 


03!^ 19™ 48.5= 


41°30'27" 


Abell 496 


1992 Sep 12 


HRI 


14488 


04'^33™38.1= 


-13°15'42" 


Virgo 


1992 Dec 17 


PSPC 


9961 


12h30™49.8= 


12°23'32" 


Centaurus 


1994 Jul 05 


PSPC 


3192 


12h48'"48.7= 


_41°18'44" 


Coma 


1991 Jun 17 


PSPC 


22108 


12h59"M6.1'= 


27°56'21" 


Abell 1795 #1 


1992 Jun 25 


HRI 


2768 


13'^48™52.7= 


26°35'27" 


Abell 1795 #2 


1993 Jan 21 


HRI 


11088 


— 


— 


Abell 1795 #3 


1994 Jun 23 


HRI 


11080 


— 


— 


Abell 2199 #1 


1991 Feb 10 


HRI 


5308 


161^28^38.5= 


39°33'03" 


Abell 2199 #2 


1994 Feb 03 


HRI 


26528 


— 


— 


Abell 2199 #3 


1994 Aug 31 


HRI 


20976 


— 


— 


Abell 2597 


1992 Jun 06 


HRI 


16228 


23'^25'"19.7= 


-12°07'27" 
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Table 4. The best-fit parameter values and 90 per cent (Ax^ = 2.71) confidence limits from the basic spectral analysis of the ASCA 
data. Temperatures (fcT), metallicities [Z), column densities (A'^h)i intrinsic column densities (AA^u)! rnass deposition rates (Mg) and 
the normalization of the intrinsically absorbed (cooler) emission components in spectral model D, were linked to take the same values 
in all four detectors. Only the normalization of the hotter isothermal emission component was allowed to vary independently for each 
detector. Temperatures are quoted in keV and metallicities as a fraction of the solar photospheric value (Anders & Grevesse 1989). 
Column densities and intrinsic column densities are in units of 10^^ atom cm"^ Mass deposition rates are in Mq yr~^ . 





Parameters 


Model A 


Model B 


Model C 


Model D 


Model E 




kT 


A ^^,+0.03 
^•^^-0.03 


4 31+004 


4.70+0-18 


6.00+0-2^ 


4 57+0-09 




Z 


47+0-02 


r, 47+0.01 


n 4fi+o-02 

^■^•^-0.02 


45+0^02 
^■^°-0.02 


0.46«-.o^ 


Abell 426 


Nn 


1.49 


-1 0-7+0.03 
-■-■•J' -0.03 


i-i4tri? 


n 00+0.12 


omtVit 




M 


— 




258tll 




303+!j:^ 




fcTa 


— 


— 


— 


9 -17+0.14 


— 




AAfH 


— 


— 


3 6+0-4 


9 c+0.5 


4+0.4 




xVdof 


2487/1659 


2443/1658 


2405/1656 


2229/1655 


2389/1656 




kT 


o 70+0. 04 


3 54+0-05 


3.86«:1^ 


oc+0.10 

'^-^'^-O.OB 


3.67+0-0? 




Z 


n t,n+o-03 


51+0-03 
^•°l-0.03 


0-50+OM 


o.5ot»:°^3 


0-50l»:»^ 


Abell 496 


Ah 


0.46 


n OO+0.05 


0A9tr2\ 


0.7410.0? 


0.44t°:g 




M 


— 


— 


ISdttl 


— 


notll 




fcTa 


— 





— 


n 47+0-08 


— 




AA^H 


— 


— 


4 0+0.6 


ii.9t?:^ 


3.6«;? 




xVdof 


1546/1283 


1400/1282 


1379/1280 


1319/1279 


1380/1280 




kT 


9 05+0-01 
^•^°-0.02 


no+O.Ol 

^•U^-0.02 


r, r|O + 0.03 

^■^•^-0.03 


9 ,q+0-02 
^■4^-0.03 


9 ns+O-02 
^•U°_0.02 




Z 


67+002 

'^•"'-0.02 


n Kc+0.02 


omt^ol 


„ 70+0.04 
u-'y-o.o3 


75+0-02 
U-'f-0.03 


Virgo 


Ah 


0.25 


n ■^5+0-04 


o.o5t°:J^ 


31+0-07 
U-d4_o.07 


n 07+0.19 




M 


— 




10.7tl:0 




5 8+1-0 




fcTa 


— 


— 




0-84in^ 


— 




ANh 


— 


— 


7^+0.28 
•^-'"-0.48 


s.oin 


-1 sS+0-53 
-l-»»_0.65 




xVdof 


2527/1153 


2502/1152 


1975/1150 


1950/1149 


2018/1150 




kT 


2 77+0.01 


a 70+0.01 
^•''J-0.02 


,„+0.07 


9^+0.08 
■^-^+-0.08 


2 92+0-03 
^•^^-0.03 




Z 


-, -,0+0.03 
^•-'-»-0.03 


1 1fi+0-02 
l-l"-0.03 


1 24+0.03 
1-^^-0.03 


f. 03+0-04 
^■^■^-0.03 


1 90+0.03 

^^•^•^-0.03 


Centaurus 


Nn 


0.81 


1 ni+0-03 

1-^1-0.03 


1 no+0-20 

1-^^-0.20 


-I 90+0. 08 

^•"'^-o.io 


1 01+0-21 




M 


— 




42.013-7 




36.3t^:* 




fcTa 


— 





— 


, 04+0.04 


— 




AA^H 


— 


— 


2 46+0-36 


n 90+0.35 

^•"=■^-0.23 


1 09+0-25 
^-"^-0.33 




xVdof 


4928/1553 


4782/1552 


2703/1549 


2584/1549 


2696/1549 
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Table 4. Spectral Results - continued 





Parameters 


Model A 


Model B 


Model C 


Model D 


Model E 




kT 


6 44+0-58 

"-^^-0.52 


°-U»_0.25 




q p,7+0.62 






Z 


p o„+0.07 


19+0-04 
^•^^-0.02 


— 


91+0-0-4 


— 


Coma 


M 


0.09 


n 24+0-05 


— 


11+0-07 


— 




fcTa 


— 


— 


— 


<0.87 


— 




AiVn 


— 


— 


— 


> 10.7 


— 




xVdof 


1696/1228 


1274/1227 


— 


1254/1224 


— 




kT 


c- 40+0-08 


r- qq + 0.10 


I- 07+0.25 


6.21IO-I 


5 61+0-17 




Z 


-Hfi+o-os 


3fi+0-03 
'^••^''-0.02 


Q 07 + 0.02 

*J-^'-0.03 


r. 07+0.03 
"-*'-0.03 


07+0.02 
^••^'-0.03 


Abell 1795 


Nn 


0.12 


n iy+0.05 


n 00+0-06 


00+0-03 


O.OOlO-O? 




M 


— 




301±76 


— 


246±^,J 




fcTa 


— 


— 


— 


1 48+0-54 

-^-4^-0.29 


— 




ANn 


— 


— 


1+0.7 
■^•4-0. 5 


5 0+2.8 


2 fi+0-7 




xVdof 


1442/1236 


1440/1235 


1419/1233 


1387/1232 


1422/1233 




kT 


4 lfi+0-05 


01 +0.06 

d-yt_o.o5 


4 -HS+O-ia 

'±.00_Q 22 


4 07+0.35 


4 09+0-08 

''•Uy.g gg 




Z 


42+0.02 
'J-4^-0.02 


43+0.02 
^•^■^-0.03 


n 43+0.03 

"-^•^-0.03 


40+0.03 

U-4U_g gg 


42+0-03 


Abell 2199 


Nn 


0.09 


n 40+0-05 


o-oot^;?g 


00+0-10 

"-UU_ggg 


o-ootl2 




M 


— 




1801?^ 




147tS 




kT2 


— 


— 


— 


1 0^+0.36 


— 




ANn 


— 


— 


4 1+0.5 
^•4-0.9 


2.St\i 


3 5+0-5 




xVDOF 


1693/1295 


1571/1294 


1550/1292 


1516/1291 


1549/1292 




kT 


o qn+0.07 


qO+0.06 


q 00+0.29 


2.31+0-^3 


3.45i;::is 




Z 


r, 09+0.03 
'-'-■'^-0.03 


qc+0.03 
U-30_o.04 


Q oc+0.05 


0-29+^:^4 


O-35+ni 


Abell 2597 


Nh 


0.25 


n qq+0-07 


o.2it°:^t 


0.65+0;5| 


0.0010-38 




M 


— 




774tjg 


— 


71 7+66 
'l'-208 




fcTa 


— 


— 


— 


6 22+1-**® 


— 




ANn 


— 


— 


r. 4+0.9 


2 4+1.7 


5 0+0-7 




xVdof 


1412/984 


1127/983 


1114/981 


1083/980 


1110/981 
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Table 5. The mass deposition rates (in Mq yr -'^ ) and intrinsic column densities (in units of lO^-' atom cm ^) determined from the fits 
to the ASCA spectra with the isothermal cooling flow models for various values of the slope parameter r}. 



Parameters 



r; = 0.75 



r] = 1.0 



r) = 1.5 



■q = 2.0 



T] = 2.5 



Abell 426 



M 
AiVn 

xVdof 



283 



+68 



74 

+0.3 

0.4 



3.2: 

2391.1/1656 



303^ 



h75 
-74 

3.4+0-4 



-0.4 

2388.9/1656 



329t®3 
2389.8/1656 



340-^ 



-,+89 
-83 

3.8+°-^ 



-0.4 

2390.5/1655 



347"* 



7+95 
-83 

3.9+°-^ 



-0.4 

2390.7/1656 



Abell 496 



M 
ANh 

xVdof 



go + 37.6 
q r, + 0.6 

•^■^-0.7 

1382.7/1280 



^^^-42 
3 6+0-6 

1379.9/1280 



136t« 

4.0t°:? 

1377.5/1280 



4.2tg;? 
1376.2/1280 



162 



-52 
-57 



4.4"* 



4+0.6 
-0.7 

1375.4/1280 



Virgo 



M 

ANn 
X^/DOF 



f-0.8 

-■-0.9 

F0.5 

-0.7 

2022.9/1150 



5.0: 

1.4^ 



c o+l.O 
^■^-0.9 
1 9+0.5 

2017.9/1150 



7.9 

2.7 



-1.2 
+0.5 



-0.6 

2004.0/1150 



)+1.3 
-1.4 

1+0.4 
-0.5 

1996.3/1150 



9.2"* 
3.1^ 



10.1 
3.3 



+ 1.4 
-1.5 
+0.4 



-0.5 

1991.5/1150 



Centaurus 



M 
ANn 

xVdof 



f2.6 

-2.0 

7+0.2 

-0.3 

2708.4/1549 



33.3; 

1.7^ 



3+2.4 
-2.8 
fO.3 
-0.3 

2696.3/1549 



36.3^ 

1.9^ 



41.7^ 



7+3.2 
-3 2 

2,3+0.3 



-0.3 

2681.6/1549 



45.3^ 



3+3.6 
-3.6 

2,6+0.3 



-0.3 

2678.2/1549 



47.7"* 



7+3.9 
-3.8 
9 7+0.3 
^■'-0.3 

2678.0/1549 



Abell 1795 



M 
ANn 

xVdof 



222tf^ 

9 q+0.7 
1422.7/1233 



2.61^:1 

1421.5/1233 



2841^^ 
1420.3/1233 



306l^«« 

3.0+oj 

1419.6/1233 



1+83 

^^-83 

1+0.7 

-0.5 

1419.2/1233 



321] 
3.2"* 



Abell 2199 



M 

ANn 
X^/DOF 



- + 17 
-40 

-0.7 

1550.1/1292 



125] 
3.2^ 



h20 
-46 
fO.5 
-0.7 

1548.7/1292 



147Z 
3.5^ 



S+22 

-54 

;j + 0.5 

-0.6 

1547.6/1292 



1687 
3.8"* 



^23 

-59 

-,+0.5 

-0.7 

1547.1/1292 



186: 
4.0"* 



7+24 

-64 

-)+0.5 

-0.7 

1546.8/1292 



197: 
4.2"* 



Abell 2597 



M 
ANn 

xVdof 



gio+59 

4 7+0.7 

1112.2/981 



71 7+66 
'-^ -208 

O.U_g 7 

1110.2/981 



7+77 
-247 

4+0.7 
-0.7 



867: 
5.4^ 
1108.4/981 



963_277 

0.0_g 7 

1107.5/981 



1029^ 



:i+91 
-299 
c 7+0.7 
^•'-0.6 

1106.9/981 
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Table 6. The results from the fits to individual SIS data sets (instrument specified in parentheses) over the 3.0 — 10.0 keV energy range, 
using spectral model A. 



Cluster 


kT 


Z 


X^/DOF 




(keV) 


(Zq) 




Abell 426 (81) 


5 77+0.46 


4fi+0-0^ 


148.1/143 


AbcU 496 (SO) 


4 14+0.38 


f, oc+O.OT 


117.7/118 


Virgo (SO) 


n 17+0.11 


c;i+0.10 


108.5/108 


Centaurus (SI) 


n lo+O.lS 

■^■-■-^-o.ie 


Q 07+0.13 
'^■^'-0.11 


146.5/142 


Coma (SO) 


o 40 + 1.04 


n 17+0.06 


146.7/149 


Abell 1795 (SO) 


'^■J^^-0.43 


n OO+0.06 


165.0/146 


Abell 2199 (SO) 


4 80+°-''2 


n 4n+o-o'' 


153.3/126 


Abell 2597 (SI) 


4 n^+0-66 


f. 04+0.13 

^■•^^-0.11 


67.6/78 
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Table 7. The best-fit parameter values and 90 per cent (Ax^ = 2.71) confidence limits from the spectral analysis with the two 
temperature model incorporating the power-law components. The normalizations of the power-law components {Ai) are quoted at an 
energy of IkeV in units of 10^'' photon keV^^cm~^s~^. The 2 — 10 keV fiuxes associated with the power-law components (-Fx.2-io) 
are in units of ergcm~^s^^ and are not corrected for absorption. The intrinsic luminosities in the 1 — 10 keV band (Lx i— lo) a^re 
corrected for absorption and are in units of ergs~^ . Ax^ values are the improvements in x^ obtained with the introduction of the 
power-law component into the best-fitting two-temperature model (model D; or in the case of the Coma Cluster the best fitting single 
temperature model B). The final x^ values and number of degrees of freedom in the fits are also listed (x^/DOF). The results for the 
Virgo and Centaurus clusters are from Allen et al. (2000). For these two systems the two temperature models incorporating variable 
element abundances (Section 4) were used in the analyses. 



Parameters 



Abell 426 



Abell 496 



Virgo 



Centaurus 



kT 
Z 

Nil 
M 
fcTa 

AA^H 

r 

Al 

^X,2-10 
-^X.1-10 

Ax^ 

xVdof 



4 4Q+0.35 
^•^'^-0.27 
2+0.08 
-0.07 



0.78^ 



0.43 



+0.19 
-0.23 



2.14^ 



1+0.45 
-0.42 
q O+0.5 



2.05 



+0.05 



-0.05 
fl007 
-1061 

X 10- 

™+0.34 
''-0.36 
171 

2058/1653 



7878 
1.80+0;- 
3.96t"it X 1044 



3.42 
0.53 
0.64 

0.50 



+0.12 
0.12 

+0.06 
0.03 
0.10 
0.16 



+0.12 
0.10 



12.3"* 



,+2.8 
-1.7 
1 44+0-44 
^•^^-0.38 
lQO+150 

10 5.58+1-42 ^ -^0-12 

3.29ti:tl X 1043 

49.7 

1269/1277 



2 01+"-°4 

^•^^-0.03 

,+0.04 

-0.04 



0.72"' 



0.30 



+0.08 
-0.08 



0.79 



+0.05 
-0.05 
4 SS+0-72 
4-»«-0.86 

-'-■^'-'-0.46 

137tlf 

s.mtlfl X 10-12 

4.22li-3^ X 1041 

90.0 

1468/1143 



O r,y + 0.26 

■^•^'-0.12 

+0.03 

0.07 



0, 



0.96 



+0.22 
-0.26 



1.46 
0.59 
0.76 
22.0 



+0.07 
0.06 

+ 1.01 
0.54 

+0.59 
0.65 

+62.1 
-16.4 

4.20tlil X 10"''' 

2.14l0:» X 1042 

38.8 

954.7/746 



Abell 1795 



Abell 2199 



Abell 2597 



kT 


5 50+0-51 
^•''"-0.88 


3.641J:- 


3.0610-.25 


o 4^+0.63 


Z 


0.36tH3 


0.441H6 


0.45«-.°^ 


20+0-06 


M 

kT2 


o.oolH? 


27+''-15 
^^•^'-0.27 


0-36tS-j^ 


f, C4+0.28 


1 OO+0.46 
J^->^°_0.45 


0-9llSi6 


0.86t»-.}« 





ANu 


7 2+5-7 
'■^-3.3 


lo.stti 


12.4t4.i 


— 


r 


o.i2+_roi 


im+Jifi 


1-80-0.66 


o q -+0.82 
■^•"^**-1.63 


Al 


5.96l««,° 


62.3^239^ 


iqo+187 


310lf,i 


-Px,2-10 


3.53+^-^J X 10-12 


7.24+3-02 ^ -^0-12 


4.21^1-25 X 10-12 


1.28tl'sl X 10-^' 


.^X.1-10 


5.96^^2^;;'^ X 1043 


3.5212-74 X 1043 


1.99J:J;ii X 1044 


8A9tfrl X 1042 


Ax2 


10.8 


7.3 


32.2 


9.4 


x2/DOF 


1376/1230 


1509/1289 


1051/978 


1265/1225 
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Table 8. The results on the individual element abundances in the central regions of the Virgo and Ccntaurus clusters using the two- 
temperature model with free-fitting element abundances (see Section 5). Error bars are the 90 per cent confidence limits (Ax = 2.71) 
on a single interesting parameter. The abundances of elements not listed in the table were linked to have the same ratios relative to their 
solar values as that of Fe. 



Parameter 


Virgo 


Ccntaurus 


kT 


r, ,0+0.04 

^■l'5_o.04 


^■09tl:% 


Nh 


n OO+0.08 


0.85+";- 


AJVh 


041+0.96 


n fi2+o-^'* 

'^•"^-0.62 


kT2 


o.86t°;°| 


1 c:t: + 0.09 

J^-^^^-o.ii 


Fe 


0.79+°:°^ 


0.93+ni 


Mg 


r, 90+0.11 

"■28-0.10 


o.sotlil 


Si 


1 r,p, + 0.09 


1.29t°il 


Na 


ll-O+sl 


— 


S 


1 nn+01'' 

^■^'-'-0.09 


1 9c:+0.24 


O 




^•^^-0.00 


X'/DOF 


468.9/357 


322.3/216 



Table 9. Results from the dcprojection analysis of the ROSAT images. Column 2 lists the binsizc (in arcsec and kpc) used in the 
analysis. Columns 3 — 5 summarize the velocity dispersions, core radii and galaxy linear mass (GLM) components (where required) used 
to parameterize the cluster mass profiles. Columns 6-8 summarize the results on the cooling time (icool ; the time for the gas in the central 
bin to cool from the ambient cluster temperature at constant pressure) , the cooling radius (r^^ol \ the radius where the cooling time of the 
cluster gas first exceeds a Hubble time) and the mass deposition rates (-Mi; the integrated mass deposition rates within the cooling radii). 
Errors on the velocity dispersions are 90 per cent confidence limits. Errors on the cooling times are the 10 and 90 percentile values from 
100 Monte Carlo simulations. The upper and lower confidence limits on the cooling radii are the points where the 10 and 90 percentiles 
exceed and become less than the Hubble time, respectively. Errors on the mass deposition rates are the 90 and 10 percentile values at 
the upper and lower limits for the cooling radius. We note that the errors on the deprojection results do not account for uncertainties in 
the total mass profiles. 







MASS PROFILES 


DEPROJECTION RESULTS 


Cluster 


binsize 


(T 


re 


GLM 


^cool 


^cool 


Ml 




(arcsec/kpc) 


(kms-1) 


(kpc) 


(MQkpc-i) 


(10° yr) 


(kpc) 


(Moyr-i) 


AbcU 426 


30/15.5 


67011° 


40 




0.741°:°^ 


188136 


5871°^ 


Abell 496 


12/10.6 


550llg 


40 


8.0 X 10"' 


0.55l°:l? 


126l®l 


llllf/ 


Virgo 


30/2.60 


45011° 


50 


8.0 X 10"' 


in+ooi 


106+11 


23.21°:^ 


Ccntaurus 


30/8.91 


43011^ 


35 


1.2 X 10" 


n 44+003 


124lt^ 


39.3l25g5 


Coma 


30/19.4 


7601^0 


300 


5.0 X 10"' 


14.2l«-| 


< 29 


< 1 


AbcU 1795 


16/26.5 


700+20 


50 


1.0 X 10" 


-, 07+0.09 

-'-••^'-0.08 


19lt?^ 


462+l0« 


Abell 2199 


16/13.6 


570ll« 


80 


1.5 X 10" 


1-oilHi 


1731.^^ 


19712? 


Abell 2597 


12/25.8 


61011° 


40 




o.86l°:J° 


162182 


4231°^ 
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Table 10. A comparison of the integrated mass deposition rates determined from the deprojection analysis, both before correction 
(Mi) and after correction for the effects of intrinsic absorption due to a uniform screen of cold gas {Mq; Section 7.1), with the values 
determined from the spectral analysis (Ms) using spectral model E (model C for the Virgo Cluster). Errors on the corrected deprojection 
results account for the statistical uncertainties in the deprojected quantities but not the errors in the the intrinsic column densities. 





INTRINSIC 


ORIGINAL 


CORRECTED 


SPECTRAL 




ABSORPTION 


DEPROJ. 


DEPROJ. 


ANALYSIS 


Cluster 


AAfH 


Mi 


Mc 


Ms 




(1021 atom cm-2) 


(Moyr-i) 


(Moyr-l) 


(Mgyr-l) 


Abell 426 


o 4+0.4 
3-4_o.4 


587t!^ 


1175t?l« 


ZO'i+ll 


Abell 496 


O R + 0.6 


llllif 


289t^^° 


110lf2 


Virgo 


o 7f.+0.28 


23.2lg;i 


47.5129-6 


io.7t;:0 


Centaurus 


1 02+0-25 
l-y^-0.33 


39.3+2Sg5 


71 q+22.6 


36.3+2-^ 


Abell 1795 


2 e+o-'^ 

^■'^-0.5 


462+108 


1212t^^^ 


246l?i 


Abell 2199 


3.5+0-5 


197tS 


599+^9 


147^20 


Abell 2597 


5.0+0;^ 


4231^1 


1159^261 


71 7+66 
'-'-'-208 



Table 11. A summary of the results on intrinsic absorption from the ASCA data. Columns 2 — 6 list the differences between the 
spectrally-determined column densities and the nominal Galactic values (in units of 102latom cm~2) for spectral models B, C, C, D and 
E. Also listed are the (90 per cent confidence) constraints on the covering fraction of the intrinsic absorber, /, determined using spectral 
model E. All results assume solar metallicity in the absorbing material. The errors on the mean values are the standard deviations for 
the sample of cooling- flow clusters (the Coma cluster is excluded). 



Cluster 


B-A 


C-A 


C'-A 


D-A 


E-A 


/ 


Abell 426 


n ia+o.03 


(5+0.4 


-0 11+0-05 

'^■11-0.04 


2.5i?,i 


4+0.4 


> 0.99 


Abell 496 


n 07+0.05 

^■■''-0.05 


4 0+0.6 
^••^-0.9 


r, 00 + 0.10 

U.3«_()06 


ii.9t?:0 


c+o.e 

d.D_() 7 


> 0.97 


Virgo 


1 0+0-0-1 

U.tU_g Q3 


0+0.3 
3-8_o.5 


n fic+o.oe 


5 0+0-6 


1 SS+0-53 
l-»«-0.65 


n cc+0.15 

0-66_o.i6 


Centaurus 


20+0-03 


n c:+0.4 
^•^-0.4 


1 r,c- + 0.06 

l-"=^_oog 


n 00+0.35 
^■^■^-0.23 


1 02+0-25 
l-y^-0.33 


> 0.98 


Coma 


15+0-05 










— 


Abell 1795 


05+0-05 


1+0.7 
•^•1-0.5 


10+0-11 


t: n + 2.8 
0-^-2 9 


2 e+o-'^ 


> 0.85 


Abell 2199 


•HI +0-05 
^■■'1-0.05 


4 1+0.5 


S9+0-09 
U-3y_o.09 


2.8+\l 


r; + 0.5 
•^■°-0.7 


> 0.96 


Abell 2597 


74+0-07 
'J- '^-0.08 


5 4+0.9 


„ 74+0.07 


2 4+1.7 
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Table 12. Column 2 lists the X-ray luminosities expected to be reprocessed into (predominantly) far infrared emission by the intrinsic 
absorbing material in the clusters. The best-fit parameters determined with spectral model E have been used in the calculations (except 
for the Virgo cluster, for which model C is preferred). Columns 3 and 4 lists the observed 60 and 100/^m IRAS fluxes within a 4 arcmin 
(radius) aperture centered on the X-ray centres (Table 3) measured using the IPAC SCANPI software and co-added IRAS scans. Error 
bars are the root-mean-square deviations in the residuals, external to the source extraction regions, after baseline subtraction. Where no 
detection was made an upper limit equal to three times the r.m.s. deviation in the residuals is given. Column 5 lists the in-scan separation 
Ai? (in arcmin) between the peak of the 100/im emission and the X-ray centre. Column 6 lists the total 1 — lOOO^m luminosities calculated 
using equation 2. For Abell 2597, the predicted reprocessed luminosity may be overestimated due to the effects of radiation damage to 
the SIS detectors (sec Section 7.2.) 
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0.63 
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Table 13. The ages of the cooling flows, in units of Gyr, determined with the three methods described in Section 8. Column 2 lists the 
measurements of Allen & Fabian (1997) from their X-ray colour deprojection study (method 1). Columns 3 and 4 lists the results obtained 
from the comparisons of the mass deposition rates determined from the spectral and image deprojection analyses, both excluding (no 
AA^h) and including (ATVjj) corrections for the effects of intrinsic absorption by a uniform screen of cold gas on the deprojection results 
(method 2). Column 5 summarizes the ages inferred from the breaks in the dcprojcctcd mass deposition profiles (method 3). 
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Table 14. The results from fits to the mass deposition profiles determined from the deprojection analysis with the broken power-law 
models. Only those data from radii interior to the 90 percentile upper limit to the cooling radii were included in the modelling, with the 
exception of the Virgo cluster where (due to uncertainties associated with the effects of the PSPC rib support structure) the analysis 
was limited to the central 78 kpc (15 arcmin) radius. Columns 2 and 3 list the slopes, »?i and »?2i of the mass deposition profiles interior 
to and external to the break radii (r]-,reak; Column 4), respectively. Column 5 lists the mean cooling time of the cluster gas at the break 
radius (the errors on the cooling time are the differences relative to the values measured at the upper and lower limits to the break radii). 
Columns 6 and 7 list the integrated mass deposition rates within the break radii, with {Mc(r < r^reak)) a-nd without (Mi(r < r^roak)) 
corrections for intrinsic absorption due to a uniform screen of cold gas. Error bars are pseudo-90 per cent confidence limits determined 
from the unweighted least-squares fits, assuming a reduced x^ value xi = l-O- 
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Figure 1. (Upper panels) The ASCA data and best-fitting models and (lower panels) residuals to the fits (in units of x) for (clockwise 
from top left) Abell 426, the Virgo Cluster, Abell 2199 and the Centaurus Cluster. For clarity, only the SI and G2 data and the results 
determined with spectral model C are shown. 
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Figure 2. The residuals from fits to the SIS spectra in the 3.0 — 10.0 keV energy range using a single-temperature emission model 
with Galactic absorption (model A). The best-fitting models have then been extrapolated to cover the full 0.6 — 10.0 keV band of 
the SIS instruments. For the cooling-flow clusters, note the clear excess in the residuals at energies between --^ 0.8 and 3.0 keV (the 
region dominated dominated by the Fe-L, Mg, Si and S line complexes) indicating the presence of gas cooler than the ambient cluster 
temperatures. The deficits in the residuals at energies below 0.8 kcV provide evidence for excess absorption. For Abell 426, the solid 
points show the results for an assumed Galactic column density of 1.0 X 10^^ atom cm~^ (Section 3.3). The dotted points show the 
results using the nominal Dickey & Lockman (1990) value of 1.49 X lO^^atom cm~^. Note also the positive residuals at high energies in 
the data for Abell 426, the Virgo Clusters and the Centaurus Cluster (see Section 4). (g) QOOO RAS, MNRAS 000, 000-000 
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Figure 3. Joint confidence contours on the abundances of (left) Mg and Fe, (centre) Si and Fe, and (right) S and Fe determined from 
the ASCA SIS spectra for the Virgo Cluster. The solar photospheric abundance scale of Anders & Grevesse (1989) is assumed. Contours 
mark the regions of 68, 90 and 99 per cent confidence (Ax^ = 2.30,4.61 and 9.21, respectively). 
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Figure 4. (a) The mass deposition rates determined from the ASCA spectra (Ms) versus the values measured from the deprojection 
analysis of the ROSAT images (Mi). No correction for the effects of intrinsic absorption on the deprojection results has been made. 
The dashed line is the best-fitting power-law model which has a slope, Q = 1.15 ± 0.33 and a normalization, P = 0.36 ± 0.65. (b) 
as for (a) but with the deprojection results corrected for the effects of intrinsic absorption by a uniform screen of cold gas with the 
best-fit parameters determined from the spectral analysis (Mq). The best-fitting power-law model has a slope, Q = 1.09 ± 0.32 and a 
normalization, P = 0.18 ± 0.37. (Errors on P and Q are la uncertainties determined by bootstrap re-sampling). 
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Figure 5. The 1 — lOOO^m luminosities determined from the IRAS 60 and lOO^m fluxes (using equation 2) plotted as a function of the 
predicted reprocessed X-ray luminosities from the cooling flows. The dotted curve is the line of equality between the values. 
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Figure 6. (Upper panels) The results from the fits with simple broken power-law models to the mass deposition profiles determined 
from the deprojection study (no correction for intrinsic absorption has been made). Only those data (marked with circles) from radii 
interior to the 90 percentile upper limit to the cooling radius were included in the fits. The vertical dashed lines mark the cooling radii 
in the clusters. The dotted lines mark the 90 percentile lower limits on the cooling radii. The lower panels show the mean cooling time 
of the cluster gas as a function of radius. 
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Figure 7. The mass deposition rates determined from the ASCA spectra (Ms) versus the values inferred from the breaks in the mass 
deposition profiles (Merk) determined from the deprojeetion study. The left and right panels show the results obtained without and with 
corrections for intrinsic absorption due to a uniform screen of cold gas, respectively. The dashed curves show the best fitting power-law 
models (see text for details). The dotted lines are the lines of equality {y = x) between the values. The Centaurus Cluster (plotted as an 
open circle) does not exhibit an obvious break in its mass deposition profile within the cooling radius and so we assume Mark = Mc or 
A^Brk = Ml from Table 10, as appropriate. 
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Figure 8. A comparison of the results on the ages of the cooUng flows from the three methods described in Section 8. The dot-dashed 
lines are the results from the X-ray colour deprojection study of Allen &: Fabian (1997; method 1). The dotted lines show the ages inferred 
from the comparison of the mass deposition rates determined from the spectral and image deprojection methods (method 2). For method 
2, the average of the absorption-corrected and uncorrected results is shown, except for Abcll 2597, for which the uncorrected result was 
unbounded and therefore, only the absorption-corrected result used. (The extent of the dotted curves mark the extrema obtained from 
the absorption-corrected and uncorrected analyses). The solid lines are the measurements based on the identifications of breaks in the 
mass deposition profiles (method 3). 
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Figure 9. The fraction of current 10^^ Mq clusters which have experienced a merger with a subcluster of mass 1 — 5 X 10^'' Mq (top to 
bottom) within a given time from the present (fi = 1.0, A = 0, 6 = 1 and a cosmic fluctuation index, n = — 1.5 are assumed.) 
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